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Abstract. We prove the complete asymptotic expansion of the integrated density of 
states of a Schrodinger operator H — — A + 6 acting in when the potential b is either 
smooth periodic, or generic quasi-periodic (finite linear combination of exponentials), 
or belongs to a wide class of almost-periodic functions. 



1. Introduction 
We consider the Schrodinger operator 
(1.1) H = -A + b 

acting in L2(M'^). The potential b = 6(x) is assumed to be real, smooth, and either peri- 
odic, or almost-periodic; in the almost-periodic case we assume that all the derivatives 
of b are almost-periodic as well. We are interested in the asymptotic behaviour of the 
(integrated) density of states A^(A) as the spectral parameter A tends to infinity. The 
density of states of H can be defined by the formula 



[1.2) N{X) = N{X-H) ■= lim ^ ' ^ ' 



L^oo {2LY 

Here, is the restriction of H to the cube [—L,L]'^ with the Dirichlet boundary 

conditions, and A^(A; A) is the counting function of the discrete spectrum of A. Later, 
we will give equivalent definitions of iV(A) which are more convenient to work with. If 
we denote by iVo(A) the density of states of the unperturbed operator Hq = —A, one 
can easily see that for positive A one has 

(1.3) iVo(A) = CA''^^ 

where 

(1-4) Orf = — and Wd = —, r-r-r 

^ ' {2'nY r(l + d/2) 
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is the volume of the unit ball in M'^. There is a long-standing conjecture that, at least 
in the case of periodic b, the density of states of H enjoys the following asymptotic 
behaviour as A — )■ oo: 

oo 

(1.5) N{X)r^X''/'(Cd + Y,e,X-^), 
meaning that for each K & N one has 

K 

(1.6) N{X) = X"'^ (C, + e,^"') + Rk{\) 

i=i 

with Rk{X) = o(A^~^). In those formulas, Cj are real numbers which depend on the 
potential h. They can be calculated relatively easily using the heat kernel invariants 
(computed in [3]); they are equal to certain integrals of the potential h and its derivatives. 
Indeed, in the paper [8], all these coefficients were computed; in particular, it turns out 
that, if d is even, then Cj vanish whenever j > d/2. 

Until recently, formula (11 .Sp was proved only in the case d = 1 in [T7] for periodic b 
and in for almost-periodic b. In the recent paper ^llj, this formula was proved in the 
case d = 2 and periodic potential. In the periodic case and d > 3, only partial results 
are known, see [2], [5], [S], [12], [IB]. In particular, in [5] it was shown that formula fll.6p 
is valid with K = 1 and R{X) = 0{X~^) with some small positive 5 when d = 3 and 
R{X) = 0(A^~ In A) when d > 3. Finally, in the multidimensional almost-periodic case, 
formula (11. 6p is known only with K = and -R(A) = 0{X~), see [IB] . 

The aim of our paper is to prove formula (11.61) with arbitrary K for all dimensions d 
and for periodic or almost-periodic potentials. In the case of periodic potential, we do 
not impose any additional assumptions (besides infinite smoothness) on it. However, if 
the potential b is almost-periodic, we need it to satisfy certain extra conditions; since 
the formulation of them requires several definitions, we will list these conditions and 
formulate our main result in the next section. 

Now we discuss the difference in the approaches of [11] and this paper. To begin with, 
let us assume that the potential b is periodic. Then we can perform the Floquet-Bloch 
decomposition (see, e.g., [13]) and express the operator if as a direct integral 

(1.7) ^ = I ^(k)c^k, 

quasi- momentum k running over 0''" - the cell of the lattice F^, dual to the lattice 
of periods F. The very first thing we need to do is to replace definition (II. 2p with 
a different one. There are two problems with definition (II. 2p . The first problem is 
that this definition is rather difficult to work with. The second problem is that this 
definition makes sense only for differential operators. And, although we are working 
with a differential operator (II. ip in the beginning, our methods require us to replace 
this operator with a pseudodifferential one, and so we need a definition which works for 
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pseudodifferential operators as well. In the periodic case, the alternative definition is 
given by the formula 

(1.8) iV(A):=-i- / N{X,H{k))dk, 

{^'^) Jot 

where N{X, H(k)) is the eigenvalue counting function of -ff(k). The first step of ob- 
taining information on the density of states is to compute the precise asymptotics of 
the eigenvalues of H(k). There are two different approaches to doing this. The first 
method, called the method of spectral projections, was developed in [10]. When using 
this method, we study, instead of H, the operator H = PjHPj, where {Pj} are spec- 
tral projections of the unperturbed operator Hq := —A. It was shown in [10] that, if we 
carefully choose these projections, then the spectra of H{]<.) and H{h) are close to each 
other. Next, we can decompose the operator H into invariant subspaces. There are two 
types of such subspaces. The first type (called stable, or non-resonant subspaces) corre- 
sponds to eigenvalues of H which are far away from other eigenvalues; in studying them 
we can use straightforward perturbation theory to compute their precise asymptotic be- 
haviour. This was done in [10], but in certain cases such computations were performed 
earlier (see, for example, [6] and [21]). The second type of subspaces (called unstable, or 
resonant) corresponds to clusters of eigenvalues of H lying close to each other. In order 
to study these eigenvalues, we have to use perturbation theory of multiple eigenvalues, 
and this theory is much more difficult and less precise than in the stable case. The 
methods of [TD] allow us to reduce the study of resonant eigenvalues to the study of a 
family of operators A + eB when £ — )■ 0. Here, A and B are finite-dimensional self-adjoint 
operators and e ~ A~^/^ is a small parameter. We are interested in the eigenvalues of 
A + eB which are perturbations of zero eigenvalues of A. Of course, we can write the for- 
mula X{A + eB) ~ ^ Aj-e-^, see [7], but the coefficients Xj will, in general, be unbounded 
functions of the quasi-momentum and, therefore, we cannot integrate these asymptotic 
expansions against dk. Paper [11] deals with this problem in the case d = 2. We study 
the operator PBP, where P is the orthogonal projection onto the kernel of A, and show 
that the cluster of eigenvalues of this operator has multiplicity at most two. This allows 
us, using the Weierstrass Preparation Theorem, to prove that the eigenvalues of A + eB 

enjoy the asymptotic formula X{A + eB) ~ ^ ^i^"'' where the coefficients 

Xj and Xj are bounded functions of the quasi-momentum k and so can be integrated 
against dk. Unfortunately, this approach does not work if d > 3, since then the cluster 
multiplicity of PBP becomes unbounded. 

The second method of obtaining asymptotic formulas for the eigenvalues of -f^(k) was 
developed in [19] and [20] and was also used in [12] . This method, which we call the gauge 
transform method, consists of constructing two pseudodifferential operators. Hi and H2- 
Here, Hi = e**ife~**, where \1/ is a bounded periodic self-adjoint pseudo-differential 
operator of order 0. Thus, the eigenvalues of Hi{k) coincide with the eigenvalues of H{k). 
The operator H2 is close to Hi in norm; also, operators H2{k) have a lot of invariant 
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subspaces. As in the previous method, these invariant subspaces can be generated by 
stable and unstable eigenvalues, and the case of the stable eigenvalues can be treated 
completely (i.e. the complete asymptotic formula for such eigenvalues can be obtained). 
The difference with the previous approach lies in the form of the restriction of the 
operator H2 to a subspace generated by a resonant eigenvalue. Let ^ be a point in 
the phase space lying in a resonant region generated by a lattice subspace QJ (see Section 
5 for the definitions and more details). Then if2, restricted to the invariant subspace 
generated by ^, has a form r^J + S'(r), where r is, essentially, the distance from ^ to 03 
and S{r) is a finite-dimensional self-adjoint operator which can grow in r, but slower than 
r^. Using the method of spectral projections, we could achieve that S{r) has a simple 
form, namely, S{r) = rA + B = r{A + eB), where e = r~^, which was the advantage of 
that method. The advantage of the method of gauge transform is that all eigenvalues of 
the reduced operator r^/ + S'(r) contribute to the integrated density of states, whereas in 
the method of spectral projections only the eigenvalues coming from the zero eigenvalues 
of A (and not even all such eigenvalues) were of interest to us. This observation makes 
the method of gauge transform much more convenient to use, despite the operator S{r) 
being more complicated than in the method of spectral projections. Indeed, the fact that 
all the eigenvalues contribute to the density of states allows us to use the residue theorem 
in order to compute the sum of contributions from all eigenvalues without computing 
the contributions from individual eigenvalues, see fllO.lip and fllO.181) . In fact, formula 
(110. ISp is the most crucial observation which has enabled us to compute the contribution 
to the density of states from the resonance regions. 

When we were working on the details of this approach, we have realized that, as a 
matter of fact, the decomposition (11.71) is not required, and all the steps can be written 
for the 'global' operators H without any references to the 'fibre' operators i?(k). This 
led us to believe that this method is likely to be applicable in a range of other settings. In 
particular, it turned out that this method works for quasi-periodic and almost-periodic 
potentials. So, let us assume that the potential b is almost-periodic. What is the analogue 
of definition (II. 8p in this case? The answer to this question can be found in [16j. One 
possible way of defining the density of states is via the von Neumann algebras; we 
discuss this approach later in our paper. However, the ultimate definition is as follows: 
let e(A;x, y) be the kernel of the spectral projection of an elliptic pseudo-differential 
operator of positive order with almost-periodic coefficients. Then, it was proved in 
Theorem 4.1 of (TBI, that the density of states of this operator satisfies (at least at its 
continuity points): 

(1.9) iV(A) = M,(e(A;x,x)), 

where M is the mean of an almost-periodic function. Our main tool in the proof will be 
formula (II. 9p . but we will use the operator-algebraic definition sometimes (for example, 
to show that density of states decreases when the operator increases, something not im- 
mediately obvious from (II. 9p ). Another useful observation which helped us in extending 
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our results to the almost-periodic case is this. Let A be an elliptic pseudo-differential op- 
erator with almost-periodic coefficients. We are usually assuming that A acts in L2(M'^). 
However, we can consider actions of A (via the same Fourier integral operator formula) in 
different vector spaces, for example in the Besicovitch space i?2(I^'^)- The space -82(1^'^) 
is the space of all formal sums 

00 

XI (x), 

where 

(1.10) e,(x) :=e^^^ 

and Yl'jLi kjP < -|-oo. It is known (see [U]) that the spectra of A acting in L2(M'^) and 
B2{M.'^) are the same, although the types of those spectra can be entirely different. It is 
very convenient, when working with the gauge transform constructions, to assume that 
all the operators involved act in B2{M.'^), although in the end we will return to operators 
acting in L2(M'^). This trick (working with operators acting in i?2(K'^)) is similar to 
working with fibre operators A(k) in the periodic case in a sense that we can freely 
consider the action of an operator on one, or finitely many, exponentials, without caring 
that these exponentials do not belong to our function space. 

It seems likely that the approach of this paper can be applied to a wider class of 
operators than (11. ip . The operators should be of the type H = Hq + b, where Hq 
has constant coefficients, and b has order smaller than Hq. We plan to consider such 
operators in a subsequent publication. 

Now we describe the structure of this paper. The proof of our main theorem consists of 
several parts, which are not always immediately related to each other; in particular, there 
is no natural order in which these parts should be presented. As a result, it is possible 
to read different sections of our paper in almost arbitrary order. The main principles we 
were following in determining the actual order of the sections were: trying to postpone 
the most difficult and technical parts of the proof for as long as possible, and trying to 
minimize the amount of references to definitions/results stated after the reference. In 
particular. Section 6 of this paper can be considered as a further general discussion of our 
approach which we have decided to postpone until the definitions and results of Sections 
2-5 have been introduced. In Section 2, we give some basic definitions, formulate the 
conditions we impose on the potential and state the main result. In Section 3, we explain 
why, instead of proving (11.61) . it would be sufficient to prove a more general asymptotic 
formula (13. ip (which includes more powers of A as well as logarithms). We also explain 
why it is enough to prove this asymptotic formula not for all large A, but only for A inside 
a fixed interval. The proof of these statements (as well as reasons why we need them) is 
similar to the corresponding section in [11]. In Section 4, we describe the definition of the 
density of states based on the operator algebraic constructions and prove several useful 
properties of N{X) which immediately follow from these constructions. In Section 5, we 
define resonance regions and prove their properties. The reader who has read several 
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of the papers [20], ^0\, P, [H], [12] may have noticed that in each of these papers the 
construction of the resonance regions is shghtly different. The reason is that each time 
we define these regions, we need to fine tune the definition taking care of the problem 
we are trying to solve. Our present paper is not an exception, and the construction 
of the resonance regions in Section 5 is different from the constructions in all papers 
mentioned above. This new construction will be extremely convenient when we are 
going to integrate the contribution from individual eigenvalues to the density of states. 
In Section 6, we describe this procedure of integrating the contribution from individual 
eigenvalues over the resonance zones in more detail. In Section 7, we introduce the 
coordinates in each resonance region (or rather we cut each resonance region into pieces 
and introduce coordinates in each piece). These coordinates are introduced so that the 
integration, described in Section 6, will be as painless as it possibly can. Each resonance 
region will have two types of coordinates. The first type is Cartesian coordinates in 2J, 
where is the quasi-lattice subspace generating the resonance region. The second set of 
coordinates is the shifted polar coordinates in 2J^. These coordinates are ideologically 
similar to the shifted polar coordinates we have introduced in [11], but the details are 
much more complicated now. Starting from Section 7, until the end of Section 10, we will 
assume that all the regions where the integration takes place are of the simplest possible 
type (the simplex case). In Sections 8 and 9, we discuss the main tool of this paper, 
the gauge transform method. A large proportion of the material contained in these 
two sections is similar to the relevant parts of [12], the only difference being definition 
(18. 3p (we need to change the norm to accommodate it to the case of almost-periodic 
coefficients) and Lemma 19.31 (this lemma was not required in [12]). In Section 10, we 
compute the contribution to the density of states from each resonance region and, first, 
reduce this contribution to the explicit integral (110. 36p and then, in Lemma 110. 4[ prove 
that this integral admits a decomposition in the powers of A and logarithms. Finally, in 
Section 11, we discuss how to reduce integration over the region of arbitrary shape to 
the simplex case. 

Acknowledgement. We are very grateful to Michael Levitin for running computer 
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to Eugene Shargorodsky for useful discussions. Thanks also go to Gerassimos Barbatis 
and Sergey Morozov who have read the preliminary version of this paper and made several 
useful comments and suggestions. The work of the first author was partially supported 
by the EPSRC grant EP/F029721/1. The second author was partially supported by the 
NSF grant DMS-0901015. 



2. Preliminaries 



Since our potential b is almost-periodic, it has the Fourier series 



(2.1) 




0G0 
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where 

(2.2) e0(x) := e^^^ 

and © is a (countable) set of frequencies. Without loss of generality we assume that 
spans R*^, and contains and is symmetric about 0; we also put 

(2.3) efc:=e + e + --- + e 

(algebraic sum taken k times) and 0oo Ufc©^ = Z{G), where for a set 5" C by 
Z{S) we denote the set of all finite linear combinations of elements in 5" with integer 
coefficients. The set ©oo is countable and non-discrete (unless the potential b is periodic). 
The first condition we impose on the potential is: 

Condition A. Suppose that 6i, . . . ,0d & ©oo- Then Z{6i, . . . , da) is discrete. 

It is easy to see that this condition can be reformulated like this: suppose, 0i, . . . ,0^ & 
©oo-Then either {0j} are linearly independent, or Yl'j=i ^j^j ~ 0' where rij G Z and not 
all rij are zeros. This reformulation shows that Condition A is generic: indeed, if we 
are choosing frequencies of b one after the other, then on each step we have to avoid 
choosing a new frequency from a countable set of hyperplanes, and this is obviously 
a generic restriction. Condition A is obviously satisfied for periodic potentials, but it 
becomes meaningful for quasi-periodic potentials (we call a function quasi-periodic, if it 
is a linear combination of finitely many exponentials). 

The rest of the conditions we have to impose describe how well we can approximate 
the potential b by means of quasi-periodic functions. In the proof we are going to work 
with quasi-periodic approximations of b, and we need these conditions to make sure that 
all estimates in the proof are uniform with respect to these approximations. 

Condition B. Let k be an arbitrary fixed natural number. Then for each sufficiently 
large real number p there is a finite set ©(A;; p) C (© fl B{^p^/^)) (where B{r) is a ball of 
radius r centered at 0) and a 'cut-off' potential 

(2-4) hk;p){^) ■■= ^ a0eg{^) 

which satisfies 

(2.5) \\b - b^k;p)\\oo < ■ 

Remark 2.1. First of all, notice that we can reformulate this condition like this: for each 
(small) a > and (small) e > there is a 'cut-off' potential b(^a,e) so that — 6(Q,;e)||oo < e 
and the frequencies of &(a;e) lie inside the ball of radius e~". However, it will be rather 
more convenient in what follows to have Condition B formulated in terms of k and p. 
This condition is obviously satisfied for quasi-periodic potentials; for periodic potentials 
it is equivalent to the infinite smoothness. For almost-periodic potentials Condition B 
does not seem to follow from the infinite smoothness of b. Note that we do not require 
the coefficients ag to be equal to the 'old' coefficients a©; indeed, sometimes one can find 
a better approximation by using procedures different than the trivial 'chopping off' of 6, 
like, for example, the Bochner-Fejer summation. 
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The next condition we need to impose is a version of the Diophantine condition on 
the frequencies of b. First, we need some definitions. We fix a natural number k (the 
choice of k will be determined later by how many terms in f ll.6p we want to obtain) and 

denote B := [&{k; p)]-^ (see for the notation) and &' := @ \ {0}. We say that QJ 

is a quasi-lattice subspace of dimension m, if is a linear span of m linear independent 
vectors 0i,...,6m with Oj G © Vj. Obviously, the zero space (which we will denote 
by X) is a quasi-lattice subspace of dimension and M'^ is a quasi-lattice subspace of 
dimension d. We denote by the collection of all quasi-lattice subspaces of dimension 
m and put V := UmV^- If ^ G M'' and 5J is a linear subspace of M*^, we denote by the 
orthogonal projection of ^ onto 5J, and put to be an orthogonal complement of 5J, 
so that ^(jj± = $, — Let QJ, ii G V. We say that these subspaces are strongly distinct, 
if neither of them is a subspace of the other one. This condition is equivalent to stating 
that if we put 211 := QJ fill, then dim 211 is strictly less than dimensions of QJ and 11. We 
put = 0(QJ, il) G [0, tt/2] to be the angle between them, i.e. the angle between 5J 2IT 
and il02n, where QJ021J is the orthogonal complement of 22J in QJ. This angle is positive 
iff and 21J are strongly distinct. We put s = s{p) = s(0) := inf sin(0(QJ, it)), where 
infimum is over all strongly distinct pairs of subspaces from V, R = R{p) := sup^^Q \0\, 
and r = r(p) := inf^^Q/ \6\. Obviously, R{p) ^ p^^'' (where the implied constant can 
depend on k and k; we say that f g if f = 0(g)). 

Condition C. For each fixed k and k the sets @{k; p) satisfying (12.41) and (12.51) can 
be chosen in such a way that for sufficiently large p we have 

(2.6) sip) > p-'l^ 
and 

(2.7) r(p) > p-^l\ 

where the implied constant (i.e. how large should p be) can depend on k and k. 

Remark 2.2. First of all, we remark that condition (12. 7p for k can be derived from 
condition (12.61) for A; + 1, but we prefer to postulate both conditions. We also note that 
Condition C is automatically satisfied for quasi-periodic potentials; for smooth periodic 
potentials Condition C is also automatically satisfied (see, for example, [10]). Finally, 
notice that condition (12.61) is equivalent to s(p) > p~"/*^ for any fixed positive a (indeed, 
this equivalence can be proved by considering sets ©([a^^A;];p) instead of @{k; p) in 
Condition B, since Condition B holds for all k). Thus, if we consider potentials of the 
form b = bper + bgua-per, whcre bper is smooth periodic and bgua-per is quasi-periodic. 
Condition C amounts to the Diophantine condition on the frequencies of bgua-per and is 
generic. 

Condition A implies the following statement. Suppose, 6i, . . . ,6i E @, I < d — 1. Let 
QJ be the span of 6i, . . . ,6i. Then each element of the set fl 03 is a linear combination 
of 6i, . . . ,01 with rational coefficients. Since the set fl 2J is finite, this implies that 
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the set Z{@ fl 5J) is discrete and is, therefore, a lattice in QJ. We denote this lattice by 
r(p;QJ). Our final condition states that this lattice cannot be too dense. 

Condition D. We can choose @{k; p) satisfying conditions B and C in such a way 
that for sufficiently large p and for each QJ G V, 23 7^ W^, we have 



(2.8) vol(2J/r(p;5J))>p- 



-i/fc 



Remark 2.3. As with Condition C, Condition D is satisfied for quasi-periodic and 
smooth periodic potentials. Also, similarly to Remark 12. 2^ condition f l2.8l) is equivalent 
to vol(Q3/r(p; 03)) > p~"/^ . Condition D is not essential for our methods and it is likely 
that this condition can be relaxed. Indeed, the only place we are using this condition is 
to get an upper bound on the number of elements in T(^), and this estimate, in turn, is 
used only to prove ( 110. ISp . However, it seems likely that there may be another way of 
proving that LHS and RHS of (110. 18p are the same. This alternative proof is more direct 
and much more difficult technically. Given that our paper is quite technically involved 
the way it is now, we have decided to present a proof which is considerably simpler, 
paying the price of assuming a slightly stronger condition on the potential. 

Remark 2.4. One final remark that concerns all conditions B-D. Given any symmetric 
set of frequencies, we can construct a real smooth almost-periodic potential h such 
that (12.11) holds, all Fourier coefficients ae are non-zero, and conditions B-D are satisfied 
(of course, the Fourier coefficients will have to converge to zero really fast). For example, 
if 5 is a limit-periodic function with Fourier coefficients going to zero exponentially, than 
all our conditions A-D are satisfied. 

Now we can formulate our main theorem. 

Theorem 2.5. Let H he an operator (11.11) with real smooth almost-periodic potential b 
satisfying Conditions A,B,C, and D. Then for each K E'N we have: 



(2.9) 



as A — )■ 00. 



Remark 2.6. Following [3], [1], and [H], it is straightforward to compute the coefficients 
Cj. For example, we have 



and 



= - W ^ ^ 



K27r) 

where M is the mean of an almost-periodic function 
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From now on, we always assume that our potential satisfies all the conditions from 
this section; we also will denote p := a/A. Given Conditions B-D, we want to introduce 
the following definition. We say that a positive function / = /(p) = /(p; k, k) satisfies 
the estimate /(p) < p°+ (resp. /(p) > p""), if for each positive e and for each k we can 
achieve /(p) < p^ (resp. /(p) > p^^) for sufficiently large p by choosing parameter k 
from Conditions B-D sufficiently large. For example, we have R{p) < p^"*", s{p) > p°~, 
'"(p) > P° ) and vol(5J/r(p; QJ)) > p°~. One can also use a standard covering argument 
to show that the number of elements in @{k; p) satisfies |0(fc;p)| < p^~^. Throughout 
the paper, we always assume that the value of k is chosen sufficiently large so that all 
inequalities of the form p""*" < p^ or p"" > p~^ we encounter in the proof are satisfied. 

Remark 2.7. As we will mention several times in this paper, the very big problem for 
both the authors and the readers is the amount of notation one has to keep in mind. 
The above definition is the first step in our desire to make as much as possible of the 
notation obsolete and eventually stop using it. 

The next statement shows a bit more how this new notation is used. 

Lemma 2.8. Suppose, 0, n^, . . . , fi^ G 0, the set {fJ-j} is linearly independent, and 
6 = J2'j=i bjf^j- Then each non-zero coefficient bj satisfies 

(2.10) P°" < |&,| <P°+. 

Proof. Let 03 G V^.i be a subspace spanned by /x^-, j = 2, . . . , d, and let e be a unit 
vector orthogonal to 03. Then the sine of the angle between 6 and 03 is {6, e)\0\~^. Thus, 
if this angle is non-zero, we have |(^,e)| > s(p)r(p) and, hence, if bi = {6,e){iJ,^,e)~^ is 
non-zero, it satisfies > r{p)s{p)R{p)~'^ > p°~. Similarly, since (Aii,e) 7^ 0, we have 
> s{p)r{p) and thus < R{p){r{p)s{p))~^ < p^"*". The proof for j 7^ 1 is 
similar. □ 

In this paper, by C or c we denote positive constants, the exact value of which can 
be different each time they occur in the text, possibly even each time they occur in the 
same formula. On the other hand, the constants which are labeled (like Ci, C3, etc) have 
their values being fixed throughout the text. Given two positive functions / and g, we 
say that f^g, or g<t^f, or g = 0{f) if the ratio j is bounded. We say f ^ g if f ^ g 
and f ^ g. 

3. Reduction to a finite interval of spectral parameter 

The main result of our paper. Theorem 12. 5^ will follow from the following theorem 
(recall that we put p := VX): 

Theorem 3.1. For each K eN we have: 

(3.1) N{p') = C,pUY. E e,,,p-^(lnp)^ + o(p-^') 

p=0 j=-d+l 
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as p ^ CO. 

Once the theorem is proved, it immediately imphes 
Corollary 3.2. For each K we have: 



K 



V i=i 

as A oo. 



Proof. First of all, we notice that [3], |1] and formula (2.9) from [16] imply that 
(3.3) 



/oo ^ 
e-*^iV(A)ciA~r(^+2)/2^g,t^^ 
j=o 

as t — )■ 0+, where qj are constants depending on the potential. Now the corollary follows 
from Theorem 13.11 and calculations similar to that of [8J. Indeed, consider the following 
integrals: 



oo 

(3.4) h{t;k,p):= I e''^ X'^ln XydX, p E Z+, k E Z; 



(3.5) 



/oo 
e-t^X-''-^\nXfdX, peZ+, k e Z. 



Elementary calculations show that 
(3.6) 



hit; k,p) = t'-' (ri-k + 1) (in^y + l^a,- (l^^)' j + fiit), 



for < 0, t > 0, 



(3.7) 

p+i p 



h{t-Xp)=t'-' (in-j +Ea;.(ln-j j +/,(t), for A: > 1, t > 0; 

(3.8) 

h{t; k,p) = t'-'^ (vi-k + ^) (in^y + (in^y j +/3(t), for any A; G Z, t > 0. 

Here, aj = aj{k,p), a'j = a'j{k,p), a" = a'-{k,p) are some constants and fj{t) = fj{t; k,p) 
are entire functions in t. Obviously, e^*^N(X)dX is an entire function in t. Compar- 
ing ([331) and (jSJ]), (jSTD, dMD, it is not difficuh to see that 

1. if (i is even then cj^p can be non-zero only if p = and j is non-positive and even; 

2. if d is odd then Cj^p can be non-zero only if p = and j is odd. 

□ 
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Thus, we can concentrate on proving Theorem 13.11 

To begin with, we choose sufficiently large po > 1 (to be fixed later on) and put 
Pn = 2p„_i = 2"po) '■= Pn', wc also define the interval /„ = [p„,4p„]. The proof of 
Theorem 13.11 will be based on the following lemma: 

Lemma 3.3. For each M G N and p & In we have: 

d-l 6M 

(3.9) N{p') = C,pUY. E e,»p-^(lnp)P + 0(p, 



n 

p=0 jr = -d+l 



Here, ej,p{n) are some real numbers depending on j, p and n (and M) satisfying 

(3.10) = 0(p(^^/3)+"). 

The constants in the 0-terms do not depend on n (but they may depend on M). The 
value of a does not depend on either n or M. 

Remark 3.4. Note that (13.91) is not a 'proper' asymptotic formula, since the coefficients 
ej^p{n) are allowed to grow with n (and, therefore, with p). 

Let us prove Theorem 13.11 assuming that we have proved Lemma 13. 3[ Let M be fixed. 
Denote 

d-l 6M 

(3.11) Nn{p'):=C,p' + J2 E e,,,{n)p-^ {In py. 

p=0 j=~d+l 

Then, whenever p E Jn '■= In-i H /„ = [p„, 2p„], we have: 

d-l 6M 

(3.12) iv„(p2)-iV„_i(p2) = 5^ hpin)p''i^^Pr, 

p=0 j=-d+l 

where 

(3.13) ij,p{'^) ■= ^j,pi'^) - ^j,pi^ - 

On the other hand, since for p E Jn have both Ni^p^) = N^^p^) + 0{p~^) and 
iV(p2) = iV„_i(p2)+0(p-*^), this implies that Y.tJoT.TLd+ihM p'K^^ Pf = 0{Pn''). 

Claim 3.5. For each j = -d + 1, . . . ,6M we have: tj^p{n) = 0{pl~^\ln p^y-^'^) . 

Proof. Put X := p^^. Then ^^Iq X]^=-d+i ^j,p(^)^''(^l)^(l^^)^ — C^iPn^^) whenever 
X G [^,Pn^]- Put y := xp„ and 



r,An) := E Q (-l)%-,s(n)(lnp„)^-^. 

s=p ^ ^ 
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Then 

d-l 6M 

(3.14) Piy)--=Y. E r,My^{\nyr = 0{l) 

p=0 j=-d+l 

whenever y G [|, 1]. Consider the following d{6M + d) functions: y^ilnyY {j = —d + 
1,...,6M, p = 0, . . . ,d — 1) and label them hi{y), ...hd{6M+d){y)- These functions are 
linearly independent on the interval [i, 1]. Therefore, there exist points yi, ...,yd{6M+d) £ 
1] such that the determinant of the matrix {hj{yi))^f^l^'^^ is non-zero. Now fl3.14p and 
the Cramer's Rule imply that for each j the values Tj^p{n) are fractions with a bounded 
expression in the numerator and a fixed non-zero number in the denominator. Therefore, 
Tj^p{n) = 0(1). This shows first that tj^d-i{n) = 0(p{~*^) and then subsequently reducing 
index p from p = d — 1 to p = we obtain tj^p{n) = 0(p:^~*^(ln PnY'^~^) as claimed. □ 

Thus, for j < M, the series Ylm=o^3<p(^) absolutely convergent; moreover, for such 
j we have: 

n oo 

e„p(n) = e,,,(0) + ^ t,,,(m) = e,,p(0) + ^ t,,,(m) + 0{pi-''{\n p^f-'-^ 

("^•l^J m=l m=l 

=: e,,p + 0{pi,''\\np^r~'-n, 

where we have denoted ej^p := ej,p(0) + Ylm=i ^i.pl'^)- 

Since ej^p{n) = 0{pn''^^^~^°') (it was one of the assumptions of lemma), we have: 

6M 

(3.16) J2 \^^A^)\Pn = 0{pt^) = 0{pn^), 

j=M 

assuming as we can without loss of generality that M is sufficiently large. Thus, when 
p E In, 'Vi'e have: 

d-l M-l 

(3.17) N{p') = Cdp' + Y, E e,,pp-^{\npr + 0{p-'^{\npr-')+0{p-^{lnpr-'). 

p=0 j=-d+l 

Since constants in O terms do not depend on n, for all p > po we have: 

d-l M-l 

N{p') = Cdp' + J2 E e,,pp-^ {In py + 0{p'f) 

p=0 j=-d+l 
p=0 j=-d+l 

Taking M = + 1, we obtain (l31|) . 

The rest of the paper is devoted to proving Lemma 13.31 We will mostly concentrate 
on obtaining formula (13. 9p . since estimate ( I3.10p will usually follow by trivial but tedious 
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arguments (like estimating coefficients in the product of several geometric series). How- 
ever, in the cases when estimating the coefficients in our infinite series would present 
difficulties, we will carry out these estimates as well. The first step of the proof is fixing 
n and fixing large k and k. The precise value of k will be chosen later in order to satisfy 
estimate (19.181) (this estimate says that the more asymptotic terms we want to have in 
(13.91) . the bigger k we need to choose; note that the choice of k does not depend on k). 
We will have several requirements on how large k should be (most of them will be of 
the form < p,^ or > p""); each time we have such an inequality, we assume that 
k is chosen sufficiently large to satisfy it. The first requirement on k we have is that 
k > M. After fixing n and k, we choose the finite set 0(fc; p„) and the approximating 
potential b(^k;p„) which satisfy all the conditions B-D. Then condition (12.51) and definition 
(11.21) imply that difference between the densities of states of operators with potentials b 
and b(^k;pn) is smaller than p~^^. Thus, from now on we will consider the operator with a 
potential &(fc;p„) and try to establish (13. 9p for this new operator. Following our policy of 
getting rid of all indexes as soon as possible (i.e. immediately after we have fixed them), 
we will denote := 0(A;;p„) and b := bi^k;pn)- This means that from now on we will 
assume that 6 is a quasi-periodic potential with its spectrum of frequencies so that 
satisfies conditions A-D with p = pn- 



4. Abstract results 

In this section,we establish several abstract results concerning density of states for 
operators with almost-periodic coefficients. In the periodic setting, these results become 
either trivial or already known, so the reader who is mostly interested in the periodic 
case, can skip this section. 

In this and further sections, we will work with pseudo-differential operators with 
almost-periodic coefficients (or symbols). These operators were studied in [15] and [T6] . 
In Section 8, we will introduce the classes of such operators. We will also see that one 
can naturally consider the action of such operators in both L2(M'^) and B2{M.'^). These 
actions have many similarities between them; in particular, the norms and (for elliptic or 
bounded operators) the spectra of operators acting in L2(M'^) and -82(1^'^) are the same, 
see [15]. As a result, often when we discuss a pseudo-differential operator with almost- 
periodic coefficients, we do not specify in which space it acts. Sometimes, however, it 
becomes important to emphasize the space where the operator acts, in which case we 
will do this. 

Following [H], we denote by 21^ a I loo factor acting in ^ := B2{W-) ® L2{W'). We 
denote by both the function e*^'' and the operator of multiplication by this function. 

is the operator of translation by ^ in L2(M'^), i.e. T^m(x) = m(x — ^). The factor 21^ is 
defined as the von Neumann algebra acting in generated by two families of operators: 



(4.1) 



{e^ ® e^, ^ G R'^} 
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and 

(4.2) {/®T^, ^ G M'^}. 

Let A = a(x, D) be a self-adjoint pseudo-differential operator with almost-periodic coeffi- 
cients such that a(x, ^) ^ 1^1"* for some m > 0. We introduce operator := a(x+y, Dy) 
acting in i^; here, x is a variable of functions in i?2(IR'^) and y is a variable of functions 
in L2(M'^). We denote by Ex{A) the spectral projection of A; by Ex{A) we denote the 
spectral projection of AK By D and T we denote the relative dimension and the relative 
trace in 21^ (see [Q]). 

If A is actually a differential operator then (see [16]) one can define the density of 
states of A (denoted by A^(A; A)) by formula (11.21) . It was also proved in [TB] that 

(4.3) iV(A; A) = T(Ea(A«)) = B{Ex{A^)^). 

Note that the relative dimension in a J/qo factor can take any non-negative value. Now 
it is natural to define the density of states for a general elliptic self-adjoint pseudo- 
dijferential operaioT with almost-periodic coefficients by (14.31) . 

By £ we denote a closed linear subspace of Jq adjoint to 21^ (see [9] for the explanation 
of the terminology) . The following lemma gives a variational description of the density 
of states. 

Lemma 4.1. 

(4.4) iV(A;A) = sup{D(i:), (A«0, 0) < A(0, 0), V0 G £}. 

Proof. By taking £ := Ex{A^)^ and using fHl3l) . we see that the LHS of (lOI) is at most 
the RHS. Suppose now that we have found a subspace £ such that D(£) > D{Ex{A^)^)- 
Then Lemma from Section VII. 37 of [9J implies that £ contains a non-zero vector 
orthogonal to Ex{A^)S^. But then {A^(f),(f)) > A(0, 0), which contradicts our assumption 
on £. This proves (14.41) . □ 

Corollary 4.2. If A>B, then N{\]A) < N{X;B). 

Corollary 4.3. Suppose, Hi and H2 are two elliptic self-adjoint pseudo- differential op- 
erators with almost-periodic coefficients such that \\Hi — H2\\ <^ Pn^'^^"^ ■ Suppose, 
N[H2]p^) satisfies asymptotic expansion (13. 9p . Then N{Hi;p'^) also satisfies (13. 9p . 

Proof. Our assumptions imply that H2 — 6 < Hi < H2 + 6, where 6 <^ pn*^^*^^"'^'*- 
Previous corollary now implies that 

(4.5) N{H2 + 5; A) < N{Hi- A) < N{H2 - 6- A). 

It remains to notice that if N{H2',p^) satisfies (13. 9p . then the difference between RHS 
and LHS of (USD is O(p-^0. □ 

Lemma 4.4. Suppose, A = a(x, D) and U = ^(x, D) are two pseudo- differential opera- 
tors with almost-periodic coefficients. Let operator A be elliptic self-adjoint and operator 
U be unitary. Then N{X; A) = N{\; U-^AU). 
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Proof. Obviously, operator [/" is unitary and {U ^AUY = (f^") ^A^UK Thus, 
(4.6) 

N{X;U-^AU) = T{Exi{U^)-^A^U^)) = T{{U^)-^Ex{A^)U^) = T{Ex{A^)) = N{X;A). 
Here, the third equahty follows, for example, from Sections 36-37, Chapter 7 of [9]. □ 

5. Resonance zones 

In this section, we define resonance regions and establish some of their properties. 
Recall the definition of the set = 0(/c;p„) as well as of the quasi-lattice subspaces 
from Section 2. As before, by 0^ we denote the algebraic sum of k copies of 0; remember 
that we consider the index k fixed. We also put 0'^ := 0^ \ {0}. For each QJ G V we put 
S"^ := {I G 23, ICI = !}• For each non-zero 6 eR'^ we put n(0) := 0\e\-\ 

Let 23 G Vjn- We say that is a flag generated by 23, if ^ is a sequence 23^ G Vj 
(j = 0, 1, . . . , m) such that 23j_i C 23^ and 23m = 23. We say that {uj}Y=i is a sequence 
generated by ^ if Uj G 23^ 23j_i and W^jW = 1 (obviously, this condition determines 
each Uj up to the multiplication by —1). We denote by 3'(23) the collection of all flags 
generated by 23. We also fix an increasing sequence of positive numbers aj {j = 1, . . . ,d) 
with «d < ^ (these numbers depend only on d) and put Lj := pn' . 

Let G 0'^. We call by resonance zone generated by 

(5.1) A(0):={|GM^ \{iMO))\<L,}. 

Suppose, G 5'(23) is a flag and {i^jj^Li is a sequence generated by ^. We define 

(5.2) A(^?):={|GM^ |(|,i.,)|<L,-}. 

If dim 23 = 1, definition (15. 2p is reduced to fl5.ip . Obviously, if C 5^2, then K{^2) C. 

Suppose, 23 G Vj. We denote 

(5.3) Hi(23) := U^^y^^)A{d). 

Note that Hi(X) = M'^ and Hi (23) = A(0) if 23 G Vi is spanned by 9. Finally, we put 

(5.4) H(23) := Hi(2J) \ (UuD«Si(il)) = Hi(23) \ (Uud^ l^sem) A(5)). 

We call S(23) the resonance region generated by 23. Very often, the region 3(j£) is called 
the non-resonance region. We, however, will omit using this terminology since we will 
treat all regions 3(23) in the same way. 

Let us establish some basic properties of resonance regions. The first set of properties 
follows immediately from the definitions. 

Lemma 5.1. (i) We have 

(5.5) U<uev 3(23) = R"". 

(a) ^ G 3i(23) Zj(f ^fjj G ^^(23), where fi(23) C 23 a certain bounded set (more 
precisely, fi(23) = 3i(23) n 23 C B{mLm) «/dim23 = m). 
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(ill) Hi(M'^) = H(M'^) IS a bounded set, H(M'^) C B{dLd); all other sets Hi(Q3) are 
unbounded. 

Now we move to slightly less obvious properties. From now on we always assume that 
Po (and thus p„) is sufficiently large. We also assume, as we always do, that the value of 
k is sufficiently large so that, for example, Ljp^^ < Lj+i. 

Lemma 5.2. Let 5J,il G V. Then (Hi(QJ) n Hi(il)) C Hi(2IT), where W := ^ + ii 
(algebraic sum). 

Proof. Assume, without loss of generality, that mi := dim 03 > dimii =: m2- If it C 53, 
then the statement of the lemma is obvious. Consider the case when 03 and il are strongly 
distinct. Suppose, ^ G (3i(Q3) fl Si(il)). Then there is a flag G 3'(Q3) such that 
^ G A(^). Let ^1 G 3^(213) be any flag such that the ffist mi elements of coincide with 
^. Let us prove that ^ G A(3^i). Let {iyj}Y=i be a sequence generated by (m = dimOU). 
Then the inclusion ^ G A(^J) implies that < Lj for j = 1, . . . ,mi. Moreover, 

a simple geometry implies |^2jj| < + |^y|)[sin(0(Q3,il))]""'^ < 2miLmiS{pn)~^ < 

-'-'mi Lmi+i- Therefore, for j > mi + 1 we have — KCan' 

Lmi+i < Lj. This shows that indeed ^ G A(^i) and, therefore, ^ G Hi(21J), which proves 
our lemma. □ 

The next statement follows immediately from Lemma [5.21 

Corollary 5.3. (i) We can re-write definition (15. 4p like this: 

(5.6) H(23):=Hi(Q3)\(Ua^ajSi(il)). 

(ii) Ifn^ii, then H(Q3) n H(il) = 0. 

(Hi) We have M'^ = U5jgvS(53) (the disjoint union). 

Lemma 5.4. Let 03 G and 53 C Q2J G Vm+i- -^ei /x be (any) unit vector from 
2U0 23. Then, for ^ G Si (03), we have ^ G Si (213) if and only if the estimate |(^, = 
I(^27-l,m)| < ^m+i /ioWs. 

Proof. In one direction the statement is obvious. Now, we assume that ^ G Si (23) fl 
Si(2IT). Let d = {2no, . . . 211^, 233} be a flag for which $, G A(^^). If 2n„ = 23 then the 
statement of the lemma is straightforward. Otherwise, we can apply the construction 
from the proof of Lemma [5.21 with il = 2IJm- This completes the proof. □ 

Lemma 5.5. We have 

(5.7) Si(Q3) n Uii25jSi(il) = Si(23) n U2B22J,dim2B=l+dim'a3l(2n). 

Proof. Indeed, obviously, the RHS of (15.71) is a subset of the LHS. On the other hand, 
suppose, il D 23 and ^ G Si (23) fl Si(il). Then ^ G A(5^) for some ^ G 3'(il). Suppose 
that 23i G 5^ is a subspace such that dim 23i = dim 23. If 23i = 23, it immediately follows 
that $, is contained in the RHS of (15. 7p . Assume that 23i ^ 23; in particular, dim 23 > 1. 
Then there exists 232 G ^ such that dim(23 + 232) = dim 23 + 1. Put 2n := 23 + 232. Since 
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£ G A(3^) C Si(5J2), by Lemma [52] we have ^ e Hi(QU), and so ^ is contained in the 
RHS of ([52]). □ 

Corollary 5.6. We can re-write (I5.4p as 

(5.8) H(QJ) := Hi(QJ) \ (UaiTD>iT,dim2n=i+dim2jSi(2IT)). 

Lemma 5.7. Let G V and G 0^. Suppose that ^ G S(QJ) and both points $, and 
^ + e are inside A(0). Then andi + E H(5J). 

Proof. If ^ 5J, then Lemma [5^ implies that ^ G Hi(Qir), where QU = span(QJ, 0), which 
contradicts our assumption ^ G H(23). 

Let us prove that ^ + G Hi(QJ). Since ^ G H(5J) C Hi(QJ), this implies that ^ G A(^) 
with ^ G 3^(2J), = {^0 = ^, . . . , QJ„ = 2J}. 

Let J be the biggest number such that 6 ^ 5Jj_i (obviously, J < m := dimQJ). We 
construct a new flag = {iio = . . . ,iim = 2J} such that 

[ X, J = 

ilj = < span(5Jj_i, 6), 0<j<J 
[^j, J > J 

We are going to prove that ^ + 6 E Let {h'j}Y=i be a sequence generated 

by di- Obviously, for j > J we have + 0,1^^) = so that |(^ + 0,i/j)\ < Lj 

if |(^, < Lj. So, assume that j < J. If j = 1, we have ui = n{Oi), so the 
assumption ^ + G A(0) implies |(^ + 0,iyi)\ < Li. Assume now that 1 < j < J. 
Then I^^J < (||^^_J + l^uJ)s(Pn)-' < 2(j - l)L,_is(p„)-i. Therefore, |(^ + 0)^,1 < 
2(j - l)%_is(p„)-i + 1^1 < 2(j - l)L,_is(p„)-i + i?(p„) < L,. Thus, |(^ + 6/,i.,)| < 
|(^ + | < Lj. This shows that, indeed, we have + 6 E A^^i) and, therefore, 
C + 6> G Hi(2J). 

Suppose now that ^ + 6 ^ "(^)- This could only happen if ^ + G Si(Q2J) for some 
2U 3 QJ. But then the previous part of the proof would imply that ^ G Si(Qir), which 
contradicts our assumption ^ G S(^). Thus, ^ + G S(5J), which finishes the proof. □ 

The next definition is almost identical to the corresponding definition from [T2] . 

Definition 5.8. Let 6, 6i, O2, . . . ,0ihe some vectors from 0'^, which are not necessarily 
distinct. 

(1) We say that two vectors ^, ry G are 6-resonant congruent if both ^ and 77 are 
inside A{0) and (^ — 77) = 10 with / G Z. In this case we write ^ O 77 mod 0. 

(2) For each ^ G M"' we denote by "^ei^) the set of all points which are 0-resonant 
congruent to ^. For 7^ we say that Tg{$,) = if ^ ^ A(0). 

(3) We say that ^ and r] are 61,62, ... , 6i-resonant congruent, if there exists a se- 
quence $,j E = 0,1, ... ,1 such that = ^5 = 'H^ ^i^d G T^^. (^j_^) for 
j = l,2,...,/. 
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(4) We say that t; G M'^ and ^ G M'^ are resonant congruent, if either ^ = r/ or ^ and 
T) are 61,62, ■■ , resonant congruent with some 61,62, ... ,61 G 0'^. The set 
of all points, resonant congruent to ^, is denoted by T(^). For points 77 G T(^) 
(note that this condition is equivalent to ^ G T(77)) we write ^7 ^ ^. 

Note that T(^) = {^} for any £ G S(X). Now Lemma [5.71 immediately implies 

Corollary 5.9. For each ^ G S(QJ) we have T(^) C S(23) and thus 

H(QJ) = U^gH(^u)T(^). 

Lemma 5.10. The diameter 0/ Y(^) is bounded above by rriLm, if ^ ^ 3(QJ), G V^- 
Proof. This follows from Lemmas 15.71 and 15.11 □ 
Lemma 5.11. For each ^ G S(5J), 2J 7^ M'^, i/ie sei T(^) is finite, and cardT(^) <^ 

Proof. This immediately follows from Lemmas 15. H 15. 7[ I5.10[ Conditions A and D and a 
standard covering argument. □ 

6. Description of the approach 

For any set C C M'^ by CP(C) we denote the orthogonal projection onto span{e^}^ge in 
B2{M.'^) and by y^(C) the same projection considered in L2(M'^), i.e. 

(6.1) y^(e) = TxcS", 

where 3^ is the Fourier transform and xe is the operator of multiplication by the charac- 
teristic function of C. Obviously, y^(C) is a well-defined (resp. non-zero) projection iff 
C is measurable (resp. has non-zero measure). We also denote S) := i?2(K'^). Let us fix 
sufficiently large n, and denote (recall that A„ = pD 

(6.2) X„ := {I G IIP G [0.7A„, 17.5A„]}. 
We also put 

(6.3) yi = yi„:=U^ex,,T(|). 

Lemma [5.101 implies that for each ^ E A we have G [0.5A„, 18A„]. In particular, we 
have 

(6.4) A n E{R'^) = 0. 
For each QJ G Vm, m < d, we put 

(6.5) yi(QJ) :=yi„nH(QJ). 
We also denote 

(6.6) A:={^^A, ||p< AJ 
and 

(6.7) A:={^^A, IIP > A J. 
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We plan to apply the gauge transform similar to the one used in [12] to the operator 
H. The details of this procedure will be explained in Sections 8 and 9; here, we just 
mention that we are going to introduce two operators: Hi and H2- The operator Hi 
is unitary equivalent to H: Hi = U~^HU, where U = e** with a bounded pseudo- 
differential operator with almost-periodic coefficients(then Lemma 14.41 implies that 
the densities of states of H and Hi are the same). Moreover, Hi = H2 + "Jl, where 
Pll < pn^'^^^^''^^ and H2 = -A + W is a self-adjoint pseudo-differential operator with 
symbol + which satisfies the following property (see section 8 for more discussion 
about pseudo-differential operators and their symbols): 

(6.8) w{e, = 0, if (^ ^ A(6I) L^eA), or {^ + 6^ A(6>) k^eA), or (0 ^ 0^). 

Now Corollary 14.31 implies that if we prove that N{p^; H2) satisfies (13.91) . then N^p^; Hi) 
(and therefore N{p'^;H)) satisfies the same asymptotic formula. This means that it is 
enough to establish the asymptotic expansion (13. 9p for the operator H2 instead of H. 
Condition (16.81) implies that for each $, G A the subspace T(T(^))i3 is an invariant 
subspace of H2 (acting, remember, in 52(M'^)); we denote its dimension by m (which is 
finite by Lemma [5. lip . We put 

(6.9) H2{rm:=H2\^^^^^. 

Note that the subspaces '?{A)S^ and '?{A)Sj are invariant as well; by H2{A) and H2{A) we 
denote the restrictions of H2 to these subspaces; we also denote by H2{A) the restriction 
of H2 to 7{A)S). Also notice that if we consider the operator H2 acting in L2(M'^), 
then 7^{A)L2{M.'^), 7^{A)L2{R'^), and T^{A)L2{M.'^) would still be invariant subspaces. 
For each ^ ^ A the operator H2{T{^)) is a finite-dimensional self-adjoint operator, 
so its spectrum is purely discrete; we denote its eigenvalues (counting multiplicities) 
by Ai(T(^)) < A2(T(^)) < •■■ < Am(f (^)) and the corresponding orthonormalized 
eigenfunctions by {/;,j_-f(4)(x)}. Next, we list all points 77 G T(^) in increasing order of 
their absolute values; thus, we have put into correspondence to each point 77 G T(^) a 
natural number t = t{r}) so that ^(77) < t{ri') if \r)\ < \r]'\. If two points 77 = (r/i, . . . ,rjd) 
and rj' = {ri[, . . . ,1]'^) have the same absolute values, we put them in the lexicographic 
order of their coordinates, i.e. we say that t{r]) < t{r}') if rji < rj'i, or rji = rj'i and 772 < 
etc. Now we define the mapping g : A ^ which puts into correspondence to each point 
rj E A the number At(r,)(T(r7)). This mapping is an injection from A onto the set of 
eigenvalues of H2, counting multiplicities (recall that we consider the operator H2 acting 
in B2(W^), so there is nothing miraculous about its spectrum consisting of eigenvalues 
and their limit points). Moreover, all eigenvalues of H2 inside the interval [0.75An, 17A„] 
have a pre-image under g. Arguments, similar to the ones used in ^[2], show that g is a 
measurable function. Similarly, we define the mapping h:A^ 52 (M'^) by the formula 
:= Then for each ^ G A the expression {'^'^)^'^J2r)&r{^) ^vi^)^viy) 
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integral kernel of the projection CP(T(^)). Therefore, we have 

(6.10) J2 Ki^)Kiy) = 1] e^We^(y). 

Another, perhaps slightly simpler way of establishing fl6.10p is just to notice that h = Fe, 
where h is a column- vector with entries {hr^}r^(^T{(,)i e is a column- vector with entries 
{e^}^g-f(^)5 and F is a unitary matrix. Then 

(6.11) 

/.^(x)/.^(y) = h(x)^h(y) = e(x)^F^Fe(y) = e(x)^e(y) = e^(x)e,(y). 

r?eT(^) r,eT(^) 

When ^ Tl, we put g{$^) := and := e^, so that now the functions g and h are 
defined on all W^. It follows from the construction that {/^-^j^gRti is an orthonormal basis 
in 52 (M'^). 

All this implies that for each A G [0.75A„, 17A„] the function 

(6.12) e(A;x,y) := (27r)-'^ / /i^(x)/^rf|, x,y G 

is the integral kernel of the spectral projection Ex{H2; B2{M.'^)) of the operator H2 in 
i?2(IR'^); here, we have denoted 

(6.13) Gx := {I G g{^) < A}. 

Notice that e(A;x, y) also gives the kernel of the spectral projection of the operator 
H2 considered in L2(M'^). Since this is the statement we will use in our proof, let us give 
a little bit more detailed proof of it. We define a mapping 



U: (27r)-'^/2 / /i^(x)/(x)rfx 



and 

(6.14) M:=A/^, 

where -H- is the equivalence relation introduced in Definition 15. 8j Lemma 15.111 and prop- 
erty (16. 4 p imply that M is measurable. It is not hard to see that ?7 is a unitary operator 
in L2(M'^) and 



: (27r)-'^/' [ h^{^)z{^)dC 
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Indeed, we have (recall the notation f l6.14p and identity f l6.10p ): 



(6.15) 



f/*f//(x) = (271^ [ [ h^{^)h^{y)f{y)dyd^ 

= {27,)'\j + [ + [) [ h{^)hMf{y)dyd^ 
J A J A J A Jw^ 

= {2n)-\f + [) [ h{^)h^)f{y)dyd^ 
JA J A JR'i 

+ (27r)-^ / / E hr,{^)hMfiy)dyd^ 

= i2n)-%[ + [) [ e^(x)^/(y)rfyrf^ 

JA JA Jr'^ 

+(27r)-^ / / E ^vi^)'^f{y)dyd^ 

Jm Jr^ ^^^(^) 

= (27r)-'^/ / e^(x)^/(y)dyrf| = /(x). 

Jr'' Jr'' 

Now we notice that for A G [0.75A„, 17A„] the function e(A;x, y) is the kernel of 
U*xG)U- Moreover, UH2{A)U* is the operator of multiplication by g acting in L2{A). 
It remains to notice that the other two restrictions of H2 satisfy UH2{A)U* < 0.75A„/ 
and UH2{A)U* > 17A„/. Now it immediately follows that for A G [A„, 16A„] e(A;x, y) 
is the kernel of the spectral projection of the operator H2 considered in L2(M'^). 

Now (14 ■3p and Theorem 4.1 from [TB] (see (11.91) ) imply the following result (note that 
since g is a measurable function, Gx is a measurable set): 

Lemma 6.1. For A G [A„, 16A„] being a continuity point of N{X; H2) we have: 

(6.16) Ar(A;iJ2) = (27r)-^vol(G'A). 

Proof. For the proof it is enough to notice that | (x) | = | (x) | = 1 for $ ^ A and 

|/i^(x)p < cardT(^) < p(f-i)°^-i+o+ 

for E Ahj (I6.10p and Lemma [5.111 and apply Lebesgue's Limit Theorem. □ 

Since points of continuity of A^(A) are dense, the asymptotic expansion proven for such 
A can be extended to all A G [A„,, 16A„] by taking the limit. Thus, our next task is to 
compute voI(Ga)- Let us put 

(6.17) i+ := G R", gi^) < p' < m 
and 

(6.18) A-:={^eR', \^\'<p'<giO}- 
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Lemma 6.2. 

(6.19) vo\{Gx) = Wdp'^ + Yo\A+ -vol A- . 

Proof. We obviously have Gx = U i+ \ i". Since A' C B{p) and i+ nB{p) = 0, 
this imphes f l6.19p . □ 

Remark 6.3. Properties of the mapping g imply that we have A~ C A. Thus, in 
order to compute A^(A), we need to analyze the behaviour of g only inside A. 

We will compute volumes of A"^ by means of integrating their characteristic functions 
in a specially chosen set of coordinates. The next section is devoted to introducing these 
coordinates. 

7. Coordinates 

In this section, we do some preparatory work before computing volA^. Namely, we 
are going to introduce a convenient set of coordinates in S(QJ). Let QJ G be fixed; 
since (1 H(]R'^) = 0, we will assume that m < d. Then, as we have seen, ^ G 3i(23) if 
and only if G fi(2J). Let {iij} be a collection of all subspaces ilj G V^+i such that 
each ilj contains 2J. Let = i-iji^) be (any) unit vector from 11^ QJ. Then it follows 
from Lemma 15.41 that for ^ G Hi(QJ), we have $, G Hi(iij) if and only if the estimate 
l(^jA*j)l = l(Caji)A*j)l ^ Lm+i holds. Thus, formula (15.81) implies that 

(7.1) H(2J) = {^e R', G n{^) k Vj K^^x, A^,.(5J))| > W}. 
The collection {yiij(QJ)} obviously coincides with 

(7.2) {nie^.), ee&j,\^}. 

The set 3(2J) is, in general, disconnected; it consists of several connected components 
which we will denote by {H(QJ)p}p^]^. Let us fix a connected component S(QJ)p. Then 

for some vectors {Aj(p)}/=i C we have 

(7.3) H(QJ)p = G G fi(5J) & Vj (C^jx, A,(P)) > 

we assume that {/ij(p)}j=i is the minimal set with this property, so that each hyperplane 
G G fi(2J) & (^2jx, Ax,(p)) = Wi}, j = l,...,Jp 

has a non-empty intersection with the boundary of H(5J)p. It is not hard to see that 
Jp > d — m. Indeed, otherwise 3(QJ)p would have non-empty intersection with Hi(QJ') 
for some 23 C 2J'. We also introduce 

(7.4) S(QJ)p := G QJ^, Vj pt^.(p)) > 0}. 

Note that our assumption that H(QJ)p is a connected component of H(23) implies that 
for any ^ G S(QJ)p and any G \ QJ we have 

(7.5) (^,0) = (^,0^,)^O. 
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We also put K := d — m — 1. 

Let us first assume that the number Jp of 'defining planes' is the minimal possible, i.e. 
Jp = K +1. We will carry on this assumption throughout most of the paper, and only in 
Section 11 we will discuss how to deal with the more general case of arbitrary S(QJ)p. If 
Jp = K + 1, then the set is linearly independent. Let a = a(p) be a unique 

point from satisfying the following conditions: (a, fXj{p)) = L^+i, j = 1, • • • , K + 1. 
Then, since the determinant of the Gram matrix of vectors is ^ p^~, we have 

|a| <S Ljn+iPn'^- We introduce the shifted cyhndrical coordinates in H(QJ)j,. These 
coordinates will be denoted by ^ = (r; Here, X = {Xi, . . . ,X„i) is an arbitrary 

set of cartesian coordinates in r2(QJ). These coordinates do not depend on the choice of 
the connected component S(2J)p. The rest of the coordinates (r, $) are shifted spherical 
coordinates in QJ-*-, centered at a. This means that 

(7.6) r(^) = l^^x - a| 
and 

(7.7) $ = n(^23x -a) G S-^jx. 

More precisely, $ e M, where M = Mp := {n^^jx - a), | G H(2J)p} C S^± is a 
X-dimensional spherical simplex with K + 1 sides. Note that 

Mp = {n(^5jx - a), ^ e E{^)p} = {n(^^x - a), Vj (|^x, /i,(p)) > L^+J 
■ = Mr]), r; := l^jx - a e Vj (r;, /x^.(p)) > 0} = S^. n H(QJ)p. 

We will denote by (i$ the spherical Lebesgue measure on Mp. For each non-zero vector 
H e ^J"*", we denote 

(7.9) Wifi):={iie^^, (r7,/x) = 0}. 

Thus, the sides of the simplex Mp are intersections of W{/lj{p)) with the sphere S<ij±. 
Each vertex v — Vt, t = I, . . . , K +1 of Mp is an intersection of Sfu± with X hyperplanes 
VF(/iy(p)), j = 1, . . . , K + 1, j t. This means that is a unit vector from QJ-*- which 
is orthogonal to {/u^(p)}, j = 1, . . . , X + 1, j 7^ this defines v up to a multiplication 
by-1. 

Lemma 7.1. Let iii and ii2 be two strongly distinct subspaces each of which is a linear 
combination of some of the vectors from {/lj{p)}. Then the angle between them is not 
smaller than s{pn) ■ In particular, all non-zero angles between two sides of any dimensions 
of Mp as well as all the distances between two vertexes Vt and Vr, t ^ t, o-f^ bounded 
below by s{pn). 

Proof. First of all, we remark that ilj are not, in general, quasi-lattice subspaces. How- 
ever, each algebraic sum QUj := QJ + ilj is a quasi-lattice subspace. Moreover, the angle 
between Wi and 22J2 is equal to the angle between ili and il2, so the first statement 
follows from Condition C. To prove the second statement, it is enough to notice that any 
non-zero angle between two sides (of arbitrary dimension) of Mp is equal to the angle 
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between two subspaces ili and il2 of the type considered in the first statement; the same 
can be said about the distance between and v,-. □ 

Lemma 7.2. Let p be fixed. Suppose, E 0^ \ and 6yj± = J2f=i ^ji^jip) ■ Then 
either all coefficients bj are non-positive, or all of them are non-negative. 

Proof. Suppose not. Then, without loss of generahty, we can assume that 6i > and 
62 < 0. Let L be a spherical interval joining vi and V2, i.e. 

(7.10) L = {ueM;, (u, il^{p)) = 0, j = 3,...,K+ 1}. 

Note that (vij^^rx) = bi{'Vi, fi^^p)) > and (v2,0ijjx) = 62(v2, /i2(p)) < 0. Therefore, 
there is a point u G L such that (u, ^5j±) = 0. This means that W{6<}ji±) has a non-empty 
intersection with Mp, which contradicts (17. 5p . □ 

Assume that the diameter of Mp is < (100(i^)~^, which we can always achieve by 
taking sufficiently large k. We put $5 := | — 0(^(jj± — a, p,g{p)), q = 1, . . . , K -\- 1. The 
geometrical meaning of these coordinates is simple: $q is the spherical distance between 
$ = n(^2j± — a) and W{fj.y{p)). The reason why we have introduced is that in these 
coordinates some important objects will be especially simple (see e.g. Lemma [775] below) 
which is very convenient for integration in Section 10. At the same time, the set of 
coordinates (r, {$g}) contains K -\-2 variables, whereas we only need K coordinates 
in QJ-*-. Thus, we have one constraint for variables $j. Namely, let {ej}, j = 1, . . . , K 
be a fixed orthonormal basis in QJ-*- chosen in such a way that the K + 1-st axis passes 
through Mp. Then we have e^- = Ylh=i^ C'jiP'i with some matrix {aji}, j,l = 1, . . . , K 
and p,i = fii{p). Therefore (recall that we denote 77 := ^^j-l ~ a), 

K+l 

(7.11) r]j = (77, e^ ) = r ^ ajg sin $g 

9=1 

and, since r^(^) = \r]\'^ = J2f=iVj^ this implies that 

(7.12) 5^(^a,,sin$,)2 = l, 

j g 

which is our constraint. 
Let us also put 

K+l 

(7.13) r/;:=^ = ^a,,sin$,. 

1^1 q=l 

Then we can write the surface element in the coordinates {rjj} as 

(7.14) d^ 



dr][ . . . dr]'j^ dr][ . . . drj'j^ 



where the denominator is bounded below by 1/2 by our choice of the basis {ej}. 
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Lemma 7.3. For each p, I we have \api\ < s{pn)~^- 

Proof. This follows from the fact that for each p \api\ is the length of the projection of Sp 
onto fii parallel to the linear space spanned by all fij, j 1. Since the absolute value of 
the sine of the angle between /i; and the hnear space spanned by all fxj, j ^ I, is at least 
s{pn), this implies that for each l,p we have \api\ < s(p„)~^, which finishes the proof. □ 

Lemma 7.4. We have maxj sin $^(77) > s{pn)d~^^'^ . 

Proof. Suppose not. Then for each / we have sin$;(r7) < s{pn)d^^^'^. Since all sin$; are 
positive, Lemma [7^ implies 

J](J]a,,sin$0' < d{ds{pnr's{pn)d-'/y = 1, 

j I 

which contradicts fl7.12p . □ 
The next lemma describes the dependence on r of all possible inner products {^,0), 

Lemma 7.5. Let $, E H(QJ)p, 2J G Vm, and 9 e 

(i) If then 0) does not depend on r. 

(a) If 6 ^ QJ and 6(^± = J^qKP'qip)' ^^^'^ 

(7.15) 6) = (X, 6>2j) + L^+i 6, + r(^) \ sin 

In the case (ii) all the coefficients hq are either non-positive or non-negative and each 
non-zero coefficient hq satisfies 

(7.16) pI < \hq\ < pI\ 
Proof. We begin by noticing that 

(7.17) {t0) = {x,e^) + {$^^,d^.), 

from which part (i) immediately follows. Recalling that ^sjj± = a + 77, we obtain 

0) = (X, 0«) + (a, 02jx) + (ry, 

= {X, 6>Qj) + V &g (a, ) + V &g (t7, l^q) 

(7.18) T Y 

= {X, e<s) + bqLra+l + T^bq Slu 

The last statement follows from Lemmas 17.21 and 12.81 The application of Lemma 17.21 
is straightforward, so let us discuss the application of Lemma 12.81 Suppose, 0ijjx = 
J2q=i^ ^gP-qip)y where 6 belongs to 0^ (but fiq{p), in general, is not in 0^). We know 
that the linear span of each fiq{p) and 5J is an element of Vm+i- Therefore, for each 
q = 1, . . . , K -\- 1, there is a vector i/g G 0^ such that p^q{p) is proportional to {iyq)<s-L, 
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= C{q){h'q)<}j±, where p° < |C(g)| < p^"*". Now we choose arbitrary hnearly 
independent vectors vk+2, ■ ■ ,i^d ^ (03 fl ©^). Then we can write = 05j± + 0<d = 

=K+2^Q^>i- Now we can apply Lemma [23] directly. □ 

8. Pseudo-differential operators 

In this and the next sections, we construct operators Hi and H2 described in Section 
6. Most of the material in these two sections is very similar to the corresponding sections 
of [l2], as are the proofs of most of the statements. Therefore, we will often omit the 
proofs, instead referring the reader to [19], [20], and [12] . 

8.1. Classes of PDO's. Before we define the pseudo-differential operators (PDO's), we 
introduce the relevant classes of symbols. 

For any / G L2(M'^) we define the Fourier transform: 

{^f m = / e-^«V(x)rfx, i G M'^. 

[2-n) 2 JiRd 

Let us now define the symbols we will consider and operators associated with them. 
Let h = 6(x, ^), X, ^ G M"', be an almost-periodic (in x) complex- valued function, i.e. 
for some countable set of frequencies (we always assume to be symmetric and to 
contain 0; starting from the middle of this section, the set will be assumed to be finite) 

(8.1) 6(x,0 = 5^&(^,Oee(x) 

6»e6) 

where 

6(0,0 :=M.(6(x,0e_,(x)) 

are Fourier coefficients of b (recall that M is the mean of an almost-periodic function). 
We always assume that (18.11) converges absolutely. Put (t) := ^1 + |t|2, Vt G R'^. We 
notice that 

(8.2) {^ + r})<2{^){ri), V^,77GM' 

We say that the symbol b belongs to the class Sq, = Sa{(3) = Sa(/3, 0), a G M, < /3 < 1, 
if for any / > and any non-negative s G Z the condition 

(8-3) \b\["J :=maxy(6>)'sup(0^""+'''^^|D|6(6l,0| < 00, |s| = si + S2 + ■ ■ ■ + s^, 

is fulfilled. The quantities (18. 3p define norms on the class Sq. Note that Sq, is an 
increasing function of a, i.e. Sq C S^ for a < 7. For later reference we write here the 
following convenient bound that follows from definition (18.31) and property (18. 2p : 
(8.4) 

^(0)'sup(O(-"+^+')^(|D|6(0,^ + ry) -D|6(0,O|) < ^ = 1^1, 



28 



L. PARNOVSKI & R. SHTERENBERG 



with a constant C depending only on a, s, and (3. For a vector t] E introduce the 
symbol 

(8.5) 6^(x,0=&(x,^ + r7),r7GM^ 

so that hrf{0^ ^) = 6(0, $, + rf) . The bound (18. 4p implies that for all |r7| < C we have 

(8.6) i&-&.i!r^<ci&iiiiN, 

uniformly in r/: \rf \ <C . 

Now we define the PDO Op (6) in the usual way: 

(8.7) Op(6)n(x) = -1- / h{-K,i)e^^-{'Jumdt 

(27r)2 J 

the integral being over W^. Under the condition 6 G Sq, the integral in the r.h.s. is clearly 
finite for any u from the Schwarz class S(]R'^). Moreover, the condition 6 G So guarantees 
the boundedness of Op (6) in L2(M'^), see Proposition 18.11 Unless otherwise stated, from 
now on S(]R'*) is taken as a natural domain for all PDO's at hand, when they act in L2. 
Applying the standard regularization procedures to definition (18. 7p (see, e.g., [15]), we 
can also consider the action of Op (6) when we apply it to an exponential Then we 
have 

(8.8) 0^{h)e^ = Y,KO.^>u+e- 

This action can be extended by linearity to all quasi-periodic functions (i.e. finite linear 
combinations of ej^). Moreover, if the order a = 0, by continuity this action can be 
extended to all of i?2(K'^); this extension has the same norm as Op(6) acting in L2 (see 
|15j). Thus, in what follows, when we speak about a pseudo-differential operator with 
almost-periodic symbol acting in B2, we mean that its domain is whole B2 (when the 
order is non-positive), or the the space of all quasi-periodic functions (for operators with 
positive order). And, when we make a statement about the norm of a pseudo-differential 
operator with almost-periodic symbol, we will not specify whether the operator acts in 
L2(M'^) or ^2(1^'^), since these norms are the same. Notice that the operator Op (6) is 
symmetric if its symbol satisfies the condition 

(8.9) hid,^) = bi-9,^ + e). 

We shall call such symbols symmetric. 

We note that in the very beginning when we consider (11. ip . our operator Op(6) is a 
multiplication by a function b (in particular, b G So). However, during modifications 
and transformations below our perturbation will eventually become a pseudo-differential 
operator. Thus, it is convenient in abstract statements to consider b a pseudo-differential 
symbol from some Sq, class. 
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8.2. Some basic results on the calculus of almost-periodic PDO's. We begin 
by listing some elementary results for almost-periodic PDO's. The proof is very similar 
(with obvious changes) to the proof of analogous statements in [19] . In what follows, if we 

need to calculate a product of two (or more) operators with some symbols bj G Sq,^. (0j) we 
will always consider that bj G S^^. (^^- 0^) where, of course, all added terms are assumed 
to have zero coefficients in front of them. 

Proposition 8.1. Suppose that \b\'^l < oo. Then Op(6) is bounded in both L2(M'^) and 
and ||0p(6)|| < 

Since Op{b)u G S(R'^) for any b e and u G S(M'^), the product Op{b)Op{g), 
b G Sq,(0i), (? G S^(02), is well defined on S(]R°'). A straightforward calculation leads to 
the following formula for the symbol bo g of the product Op(6) Op{g): 

(6o^)(x,|) = Yl &(^,l + <^)M'/>,l)e^^"+'^)^ 

ee&i,(i)e&2 

and hence 

(8.10) (T^)(x, J2 ^ + '^)^(^' 0, X e 01 + 02, ^ G M". 
We have 



Proposition 8.2. Let b G Sa(0i), g G S^(02). Then bo g e Sq+^(0i + 02) and 

1(7) 

li+(|Q|+s)/3,s' 



\bo9\\T'^<C\b\\^\g\\^^ 



with a constant C depending only on /, a, s. 

We are also interested in the estimates for symbols of commutators. For PDO's 
A,^i, / = 1,2, . . . , A^, denote 

ad(v4; ^1, ^2, . • . , ^7v) = ^ [ad(v4; ^1, ^2, • • • , ^7V-i), ^7v] , 

ad(v4; = i[A, % ad^(A; ^) = ad(A; ^, ^, . . . , ^), ad°(A; ^) = A. 

For the sake of convenience we use the notation ad(a; ijji, ip2i ■ ■ 1 "^n) and ad^(a, ip) for 
the symbols of multiple commutators. It follows from fl8.10p that the Fourier coefficients 
of the symbol ad(6, g) are given by 

(8.11) ad(M)(x,0 = ^ E [KO,i + ct^)g{4>,i)-KO,i)g{4^,i + e)\, i e 

0+<l>=x 

Proposition 8.3. Let b G Sq(0) and gj G S^^. (0j), j = 1, 2, . . . , A^. Then ad(6; gi, . . . , g^) 
S^(0 + ^^.0j) with 

N 

7 = a + ^(7^. -1), 
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N 



1(7.) 



(8.12) I ad{b■,g^, . . . ,5iv)|!J < C|6|(^)^^n M^Z^n-.^i^ 

where C and p depend on I, s, N, a and 7^. 

8.3. Pcirtition of the perturbation. Prom now on we fix ^ : < /3 < cti, and put 
Q :— & which is finite. The symbols we are going to construct will depend on p„; this 
dependence will usually be omitted from the notation. 
Let L e C°^(M) be a non- negative function such that 



.13) 



< t < 1, l{z) 



0, z>^. 



For & &,0 ^ 0, define the following C°°-cut-off functions 

\^ + 9/2\-3p, 



(8.14) 



and 



(8.15) 



lOpn 

\^ + e/2\-3pn 

lOpn 

3pn-\^ + e/2\ 
{e,^ + e/2)\ \ 



Note that ee + e^+i^ ^1. The function is supported on the set |^ + 0/2| > llp„/2, 
and ig is supported on the set \^ + 6/2\ < p„/2. The function eg is supported in the 
shell p„/4 < \^ + 6/2\ < 23p„/4. Using the notation ig for any of the functions ig or 



Eg, we point out that 



^.16) 



'eeiO = e_,(| + d), 4(|) = £_,(| + d), 

M^)^^M^ + 0), Ce(O = C-e(C + 0)- 

Note that the above functions satisfy the estimates 

f|D|e,(OI + |D|4(€)|«Pn'^ 
\|D|(^,(0| + |D|C«(OI«Pn' 



(8.17) 
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Using the above cut-off functions, for any symbol b G Sq,(/3) we introduce five new 
symbols b^^, b°, b^^, 6-^^, b^ in the following way: 



(8.18) 






(8.19) 


6^^(x,^;p„) = 




(8.20) 


6^(x,^;p,) = 




(8.21) 


&''(x,C;Pn) 


6»e0' 


(8.22) 


b°{x,$,;pn 





The superscripts here are chosen to mean correspondingly: 'large energy', 'non- resonant', 
'resonant', 'small energy' and 0-th Fourier coefficient. The corresponding operators are 
denoted by 

= Op(6^^), = Op(6^^), 
5^ = Op(6^), = Op(6''), B" = Op(6°). 
By definitions (l87[3i) . (18141) and (18351) 

The role of each of these operator is easy to explain. Note that on the support of the 
functions b^^{6, ■ ; pn) and b^{6, ■ ; pn) we have 



(8.23) 



1^1 < p:\ \pn <\i + e/2\ < |p„, lp„ - ip°+ < 1^1 < ^p„ + Ipr. 



4' 



4 



On the support of b {6, ■ ; p„) we have 



(8.24) 



< -Pn, 1^1 < -Pn + T^Pn- 



On the support of b (0, ■ ; pn) we have 



(8.25) 



11 1^1 11 1 n+ 

> yPn, III > —Pn - -Pn- 



The introduced symbols play a central role in the proof of Lemma 13.31 As we have seen 
in Section 6, due to (I8.24p and (18.251) the symbols b^^ and b'~'^ make only a negligible 
contribution to the spectrum of the operator H near the point A = p^, p & In- The only 
significant components of b are the symbols b^"^, b^ and b°. The symbol b° will remain 
as it is, and the symbol 6"^^ will be transformed in the next Section to another symbol, 
independent of x. 
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We will often combine B^, B^^ and B^^: for instance B^'^^ = B^ + B^^, B^'^^^^^ = 
j^%LE _^ j^sE. _ A similar convention applies to the symbols. Under the condition b G Sq(/3) 
the above symbols belong to the same class Sa(/3) and the following bounds hold: 

(8.26) + + |6^^|JJ + + « 

Indeed, let us check this for the symbol b^^, for instance. According to (18.231) and (18.171) . 
on the support of the function b^^{6, ■ ; p„) we have 

|DV.(^)|«P^^I^I«(0"'^'^ 

+ + |D%(OI « p;l^l « (0"'^'^- 

This immediately leads to the bound of the form (I8.26P for the symbol b^^. 

The introduced operations also preserve symmetry. Indeed, let us calculate using 



b^-e, ^ + 0) = bi-0, ^ + 0)C-ei^ + e)e^e{i + 6) 

Therefore, by ( 18. 9p the operator B"^ is symmetric if so is Op (6). The proof is similar for 
the rest of the operators introduced above. 

Let us list some other elementary properties of the introduced operators. In the lemma 
below we use the projection T(C), C C M whose definition was given in Section 6. 

Lemma 8.4. Let b G Sa{P) with some a G M. Then the following hold: 

(i) The operator B^^ is hounded and 

Moreover, 

(J - T(5(2p„/3))5^^ = B^^{I - T(S(2p„/3)) = 0. 

(ii) The operator B^ satisfies the following relations 

y(B(p„/8))i?^ = S^y(S(p„/8)) 

(8.27) = (/ - 7{B{Qpn))B'' = B''{I- T(fi(6p„))) = 0, 

and similar relations hold for the operator i?-^^ as well. 
Moreover, for any 7 G M one has b^^, b^ G and 

(8.28) l&^^lSJ+|fe^lg«p(^^"-^^l&i;?, 

for all I and s, with an implied constant independent ofb and n > 1. In particular, 
the operators i?-^^, B^ are bounded and 

115^^11 + 115^11 «pf|6|(^). 
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(iii) 

Proof. Proof of (i). It follows from (EJ]) and IK2^ that 

(8.29) 5^sup \b'\0,tpn)\ «pr^^("'°)^sup(0-^1&(^,^;Pn)| = |&lSpr^^^"'°^- 

By Proposition 18.11 this implies the sought bound for the norm 

In view of f l8.24p . the second part of statement (i) follows from fl8.2ip . 

Proof of (ii). Relations (lOTD follow from definitions (ESQ]) and (ISTTOD in view of (K23h . 

Furthermore, by (Km and (Km . 

J2{0y sup (^)(-^+^)^|D|6^=^(0, e p„)| « p^("-^) sup (O(-"+^)^|D|6^^(0, ^; p^j 

e ^ e ^ 

This means that b^^ G for any 7 G M and fl8.28p holds for b^^. The bound for the 
norm follow from (18.281) with 7 = 0, and Proposition 18. II The proof for is analogous. 
Proof of (iii) is similar to (i). The required result follows from (I8.25p . 

□ 



9. Gauge transform and the symbol of the resulting operator 

9.1. Preparation. Our strategy will be to find a unitary operator which reduces H = 
Hq + Op(6), Ho := —A, to another PDO, whose symbol, essentially, depends only on ^ 
(notice that now we have started to distinguish between the potential b and the operator 
of multiplication by it Op (6)). More precisely, we want to find operators Hi and H2 with 
the properties discussed in Section 6. The unitary operator will be constructed in the 
form U = e** with a suitable bounded self-adjoint quasi-periodic PDO This is why 
we sometimes call it a 'gauge transform'. It is useful to consider e** as an element of the 
group 

U{t) = exp{im}, Vt G M. 

We assume that the operator ad(ifo) ^) is bounded, so that U{t)D(^HQ) = D(^Hq). 
This assumption will be justified later on. Let us express the operator 

At := U{~t)HU{t) 

via its (weak) derivative with respect to t: 

At = H+ [ U{-t')8id{H;^)U{t')dt'. 
Jo 
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By induction it is easy to show that 

k 

J 



(9.1) A, = H + J2\ ad^ (^; ^) + 4+1' 

Jo Jo Jo 



'0 ^0 JO 

The operator \t' is sought in the form 

k 



(9.2) ^ = E'^^' ^^• = 0P(^^- 



with symbols ipj from some suitable class 8^^. to be specified later on. Substitute this 
formula in (19. ip and rewrite, regrouping the terms: 

j=l i=i A:i+fe2+---+fcj=; 

k 

(9.3) 4+1 ^=E^ E ad(//;vl/,^,vl/,,,...,vi/,j. 

i=l A;i+A:2H hA:j>fc+l 

Changing this expression yet again produces 
A, = H^ + Op(6) + ad(i^o; + E " E E ^^(^O' '^'^^ 



1=1 j=2 l=j fci+fc2H \-kj=l 



fe ^ k 



+ E 7f E E ^d(op(6); VP,,, VP.,, . . . , VP.J + 4^)^ + 4^_;^. 

j=l 1=3 ki+k2+---+kj=l 

Next, we switch the summation signs and decrease / by one in the second summation: 

k it I 



A^ = Ho + Op(6) + ad(i^o; ^0 + E E ^ E ^^(^o; ^a.^, • • • , ^fc,) 

1=1 1=2 j=2 ^' ki+k2 + ---+kj=l 

+ E E 7^ E ^d(op(6); VI/,,, VP,,, ... , vp,j + + 4^)^. 



1=2 j=l fci+feH hfcj=i-l 
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Now we introduce the notation 

i?i := Op(6), 

i-i 

(9.4) Bi:=J2- E ad(0p(6); vp^,, v]/,^, . . . , v|/,J, / > 2, 

j=l ki+k2+---+kj=l-l 
' 1 

(9.5) Tr.= J2- Yl ^^(^o; ^k„^k„ . . . , ^fc,), / > 2. 

i=2 ki+k2+---+kj=l 

We emphasise that the operators Bi and T; depend only on \l/2, . . . , "^i-i- Let us 
make one more rearrangement: 

k k k 

A^ = Ho + Op(6) + J2 ad(i^o, ^0 + 5/ + + ^^+1, 

(9-6) R-^^i = Bui + Rtl + R?l 

Now we can specify our algorithm for finding "^/s. The symbols tpj will be found from 
the following system of commutator equations: 

(9.7) ad(i/o; ^i) + = 0, 

(9.8) ad(/7o; ^0 + + ^z'"'' = 0, / > 2, 
and hence 

yk = j:LBi + ELTi- 

Below we denote by the symbol of the PDO F^. Recall that by Lemma I8.4( ii). 
the operators Bf'^,Tl'^ are bounded, and therefore, in view of (19.71) . (19.81) . so is the 
commutator ad(ifo;^)- This justifies the assumption made in the beginning of the 
formal calculations in this section. 

9.2. Commutator equations. Put 

~ td\ (^\(\^ , /3i2 iti2A-i (^e{$,)'^e{$,) 

Xe[i) ■■= e9[i)ve[i){\i + 0\ - \i\ ) g, 

when 7^ 0, and Xo(^) = 0. We have 

Lemma 9.1. Let A = Op(a) be a symmetric PDO with a G S^^. Then the PDO \1/ with 
the Fourier coefficients of the symbol ipl'x,^) given by 

(9.10) ^{e,^) = td{e,^)xe{^) 

solves the equation 

(9.11) ad(i/o; ^) + Op(a^^) = 0. 



(9.9) 
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Moreover, the operator \1/ is bounded and self-adjoint, its symbol belongs to with 
any 7 G M and the following bound holds: 

(9.12) l^lSJ « P^'-'-MPn)-da\tXs « P^n'^-'-'^^'^-Mtis- 

Using Propositions I8.1f8.2f8.3[ Lemma I9.H and repeating arguments from the proof 
of Lemma 4.2 from [T2j, we obtain the following estimates for the symbols introduced 
above: 

Lemma 9.2. Let b G So(/3) be a symmetric symbol. Then ipj, bj, tj G S^(/3) for any 
7 G M and 

(9.13) < c.pT-'-'MPnV (\h\t)' , J > 1; 

(9.14) |6,|SJ + |t,|SJ < C.p^^^'-^-'HpnV^' [\b\fX)' , J > 2. 

Here Cj, Ij, Sj depend only on j,l,s and 7. Moreover, assuming po is large enough 
( depending on I, s, 7, b and k) we get 

(9.15) mts^'^Pn'^'-'MPnr'lbC, 

(0) ^ oim(o). 



(9.16) \ydu<mu-, 



fc+1 



(9.17) l|i?fc+ill«Pn''^r(p„)-^(|6|!°|^,,^^^^ 
Now, we take 

(9.18) k> {M +{d-2))/P 

and assume that k is large enough so that r(p„)~^ ^ p°+ ^ p^. Then 

and we can disregard -R^^_]^ due to Corollary 14.31 More precisely, let W = Wj, be the 
operator with symbol 

(9.19) ^^(x,|):=y^(x,0-yf(x,0, i.e. t&^(0, = yfc(e, 0(1 " e«(^)^«(0). 

We put Hi := A\ and H2 := —A + W . Then — H2\\ ^ Pn*"^^*"^ and, moreover, 
the symbol w satisfies condition (16. 8p . This means that all the constructions of Section 
6 are valid, and all we need to do is to compute voI(Ga). 
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9.3. Computing the symbol of the operator after gauge transform. The follow- 
ing lemma provides us with more explicit form of the symbol yj^. 

Lemma 9.3. We have y^{0,^) = for 6 0^. Otherwise, 

k-l s 

(9.20) 



where the second sums are taken over all Oj G ©, 6'-, 0' G ©s+i and 



p 



(9.21) C.(0,O = E E C'(^H^)n^«;'(^ + '^>«"(^ + '^p- 

p=i e'^',4,'Je&s+i (i<j<p) i=i 

i/ere ci^\6) depend on s, p and a// vectors 6,6j^6'-^(j}'-^6"-^(j)"-. At the same time, 

coefficients C^f\6) can he hounded uniformly hy a constant which depends on s only. 
We apply the convention that 0/0 = 0. 

Proof. We will prove the lemma by induction. Namely, let £ > 2. We claim that: 

1) For any m = 1, — 1, ipm{G, ^) = for ^ ©„. Otherwise, 

m 

(9.22) 4^U0, ^) = Y. C'mi^^ n ^(^^•)*^^ + ^^■)' 

i=i 

where the sum is taken over all 6j G ©, 0'^, 4>'j G ©m and C'^{0, ^) admit representation 
similar to flOTD . 

2) For any s = 1, 1 and any ki, . . . ,kp {p > 1) such that ki + ■ ■ ■ + kp = s, 
ad(0p(6)i^. . . , ^fcj(6>,^) = for 6> ^ ©,+1. Otherwise, 

s 

(9.23) ad(0p(6)^ . . . , v&.J (0, ^) = ^ c;(0, mOs+i) J] K^.Ox^?; + <f>'j), 

i=i 

where the sum is taken over all dj G ©, O'j, <p'j G @s+i and C'^{6, ^) admit representation 
similar to flOTD . 

3) For any s = 2, . . . ,i and any ki, . . . ,kp {p > 2) such that ki + ■ ■ ■ + kp = s, 
ad{Ho; • , ^k,){0, C) = for 6> ^ ©,. Otherwise, 

s-l 

(9.24) ad{Ho; CT- • , ^fcj(^, ^) = (0, ^)&(0.) J] + 0;), 

where the sum is taken over all 6j G ©, O'j, <j)'j G ©s and C'"{6,$,) admit representation 
similar to (ICTD . 
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For £ = 2 statements l)-3) can be easily checked. Indeed, 

(9.25) MO,^)=tb{e)xe{^), 

(9.26) ad{6^;^^){e,0= + x) " KxW)x<^(l)) , 

X+4>=9 

(9.27) ad(i^,vl/O(0,O = - {b''''ix)Km4^ + x)-b''''{x)b{m4^))- 

x+4>=0 

Now, we complete the induction in several steps. 

Step 1. First of all, notice that due to f l9.4p . fl9.5p . for any m = 2, . . . ,i symbol of 
Bm admits representation of the form 09.231) with s = m — 1, and symbol of admits 
representation of the form fl9.24p with s = m. Then it follows from Lemma [9TT] and (19. 8p 
that admits representation of the form (19.220 . 

Step 2. Proof of flOSD with s = i. Let h + ■■■ + kp = i. lfp>2 then 

ad{Op{by, ^fc,, . . . , = ad{ad{Op{by, ^fc,, . . . , ^fc^_J; ^^J. 

Since fci + ■ ■ ■ + fcp-i < i — 1 and fcp < £ — 1 we can apply (I9.22p and fl9.23p . Combined 
with (18. lip it gives representation of the form (19.230 . If p = 1 then ad{Op{b)] \l/£) satisfies 
(19:231) because of IKTT\f and step 1. 

Step 3. Proof of flOD with s = i+l. Let h + ■■■ + kp = i + 1, p > 2. If p > 3 then 
(cf. step 2) 

ad{Ho; ^'fci, . . . , ^k,) = ad{ad{Ho; ^'fc^, . . . , ^k,.,); ^fcj- 

Since fci H h kp_i <£, p - 1 > 2 and kp < i - 1 we can apply fl^:^ and dlSlD- 

Together with (18.110 it gives representation of the form (I9.24p . If p = 2 then (see (19. 8p ) 

ad{Ho; ^fc,, ^fcj = adiad{Ho; ^^J; ^^J = -ad{B^^ + T^^^; ^^J. 

Since ki < i and /c2 < £, the representation of the form (19.240 follows from (18. lip and 
step 1. (Formally exceptional case ki = 1, k2 = i can be treated separately in the same 
way using (19.70 instead of (19.80 .) 
Induction is complete. 

Now, ([123]), and ([131), (S, <B> P^ove the lemma. □ 

10. Contribution from various resonance regions 

10.1. Summing the contributions from individual eigenvalues. Let us fix a sub- 
space QJ e Vm, m < d, and a component Hp of the resonance region H(5J). Our aim is to 
compute the contribution to the density of states from each component Hp. This means 
that we define A+(Hp) := fl H(QJ)p and A~(Hp) := A~ fl H(5J)p and try to compute 

(10.1) voli+(Hp) - voli"(Hp). 
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Since formulas f l6.19p and (16. 4p obviously imply that 

d-i 

(10.2) voKGa) = Wdp'^ + Y, Yl $^Mi+(Sp) - voli-(Hp)), 

m=0 2JGV„i p 

if we manage to compute (110.11) (or at least prove that this expression admits a complete 
asymptotic expansion in p), Lemma [331 would be proved. Thus, we fix QJ and, moreover, 
we fix a component S(QJ)p of the resonance region. Recall that K = d — m — 1. 
Note that if ^ G 3p, then we also have that T(^) C Sp. We denote 

as an operator acting in i^^ := 7{Y{^))S} (recall that S)^ is an invariant subspace of H2 
acting in i?2(K'^)). Suppose now that two points ^ and rj have the same coordinates X and 
$ and different coordinates r. Then ^ G 3p implies 77 G 3p and T{r]) = T(^) + {t] — 
This shows that two spaces Sj^ and i^^ have the same dimension and, moreover, there is 
a natural isometry F^ rj : — > -^rj given hj F : e^^(^r)-^), £ f (^)- This isometry 

allows us to 'compare' operators acting in S^^ and Sjr,- Thus, abusing slightly our notation, 
we can assume that H2{$,) and H2{r]) act in the same (finite dimensional) Hilbert space 
S^{X,^). We will fix the values (X, $) and study how these operators depend on r. 
Thus, we denote by H2{r) = H2{r; X,^) the operator -f^2(^) with ^ = (X, r, $), acting 
inSj{X, $). 

As we have seen from the previous sections, the symbol of the operator H2 satisfies 

(10.3) /i2(x, ^) = + u;^(x, = + 2r (a, 11(77)) + + W',{^, + 1^1', 

where the Fourier coefficients of Wj^ satisfy (19.161) . (19.191) . (I9.20p . (I9.2ip and we denote, as 
usual, 77 = ^jrjx — a. This immediately implies that the operator H2{r) is monotonically 
increasing in r; in particular, all its eigenvalues Xj{H2{r)) are increasing in r. Thus, 
the function g{^) (defined in Section 6) is an increasing function of r(^) if we fix other 
coordinates of ^, so the equation 

(10.4) g{0 = P' 

has a unique solution if we fix the values (X, $); we denote the r-coordinate of this 
solution by r = r(p) = r(p; X, $), so that 

(10.5) g{^iX,T,^) = p'. 

By To = To(p) = To(p; X, $) we denote the value of r for the unperturbed operator, i.e. 
To is a unique solution of the equation 

(10.6) |^(X,ro,l>)|=p. 

Obviously, we can write down a precise analytic expression for tq (and we have done 
this in [11] in the two-dimensional case) and show that it allows an expansion in powers 
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of p and In p, but we will not need it. The definition of the sets implies that the 
intersection 

(10.7) i+ n {^(X, r, $), re M+} 

consists of points with r-coordinate belonging to the interval [ro(p),T(p)] (where we 
assume the interval to be empty if tq > r). Similarly, the intersection 

(10.8) A- n {^(X, r, $), re M+} 

consists of points with r-coordinate belonging to the interval [r(p), ro(p)]. Therefore, 

(10.9) i+(Hp) = {^ = ^(X,r,<l), XGfi(5J),l'eMp, r G [ro(p; X, $), r(p; X, $)]} 
and 

(10.10) i-(Hp) = {^ = C(X,r,|.), Xen{^),~^eM„ r G [r(p;X,l>),ro(p;X,$)]}. 
This implies that 

r f ~ r^P'^^'^) 

voli+(Hp) - voli"(Hp) = / dX d^ r^dr 

= {K + l)-i f d^ f dX{T{p- X, <l)^+i - ro(p; X, <l)^+i). 

Obviously, it is enough to compute the part of fllO.lip containing r, since the second 
part (containing tq) can be computed analogously. We start by considering 

(10.12) [ T{p;X,^f+^dX. 

First of all, we notice that if ^, 77 G S(QJ) are equivalent points then, according to Lemma 
15.71 all vectors Oj from the definition [53] of equivalence belong to QJ. This naturally leads 
to the definition of equivalence for projections and t;^;. Namely, we say that two points 
u and n from f2(QJ) are QJ-equivalent (and write u /x) if u and fx are equivalent in 
the sense of Definition 15.81 with additional requirement that all 6j G 2J. Then ^ rj 
implies rj^. For u G ^(QJ) we denote by Tqj(^') the class of equivalence of u 

generated by ^^j- Then Tqj(^2j) is a projection of T(^) to QJ and is, therefore, finite. 

Denote by M^j the quotient space Msu := fi(5J)/ Since Tfa(i/) is a finite set for 
each e ^^(2J), there is a natural measure on Maj generated by the Lebesgue measure 
on fi(5J). Therefore, we can re-write (110. 12p as 



(10.13) / r(p;X,<|.)^+W 

and try to compute 



(10.14) r(p;X,$ 
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Let us denote hj S = S{r) the operator with symbol 2r(a, n(?7)) + |ap + w^(x, + |Xp 
acting in S){X, $), so that H2{r) = r'^I + S{r). 

Remark 10.1. We always assume that ^ E A, so that 0.7 pn < |^| < 5p„ and all 
functions eei^ + 0) from flOJOD - flOTD are equal to 1. Note that if G 0^, G 
and ^ QJ, then (see Lemma [731 and flS.lSp ) (P0{^ + <p) = 1- This means that all cut-off 
functions from fl9.19p - fl9.2ip are equal to 1 unless G 5J. If, on the other hand, G QJ, 
then <^6»(^ + 0) depends only on the projection and thus is a function only of the 
coordinates X. Thus, equations f l9.19p -f l9.2ip show that H2{r) depends on r analytically, 
so we can and will consider the family H2{z) with complex values of the parameter z. 

Formulas fl919D . flOO]) imply 

(10.15) ll^(r)||«p;.+"''+°+, \\S'ir)\\<^p:^^'\ 
and 

(10.16) II^^WII«P;', l>2. 

Let 7 : — p| = Pn/8} be a circle in the complex plane going in the positive direction. 
Then for p G /„ all r(p; X, $) lie inside 7. It is not hard to see that estimates fllO.lSp 
and fll0.16p hold inside and on 7 (indeed, formulas fl9.19p - fl9.2ip give matrix elements of 
S{z) in an orthonormal basis even for complex z). 

A version of the Jacobi's formula states that for any differentiable invertible matrix- 
valued function F{z) we have 

tT[F'{z)F-\z)] = {det[F{z)]y{det[F{z)])-^ 

(it can be proved, for example, using the expansion of the determinant along rows and 
the induction in the size of F) Let #(5*, 7) be the total number of zeros (counting 
multiplicity) of det[S'(2;) + z"^! — p^I] inside 7. We have 

#(5,7) = -i^ /(det[5(z) + zH - p2/])'(det[5(z) + zH - pH])-^dz 
I'm 

(10.17) 
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= ^ tr[(2^/ + S\z)){S{z) + zH - pH)-'\dz = tr[(l + 0(p^-i+°+))/], 

where / = Since cardTsjj(^') <^ pn ^^"''^^'-^^^ by Lemma [5.111 and daa < 1, we 

conclude that there are precisely cardTfj7(i/) zeros (counting multiplicities) of det[S'(2;) + 
z'^I—p^I] inside 7, and thus the points r(p; X, $) are the only zeros of det[S'(2;)-|-2;^/— p^J] 
inside 7. 



42 



L. PARNOVSKI & R. SHTERENBERG 



Then we have by the residue theorem: 



XeTa,(i/) 
1 

1 



(10.18) 



— (I z''+\det[S{z) + z^I - p^I]y{det[S{z) + z^I - p'l]Y^dz 

— ItT[z''+\2zI + S'{z)){S{z) + zH - p^iy']dz 

— / tr[(2z^+2j _^ ^i^+i5'(^))(^2 _ p2yi J^^_iys\z){z' - p')-']dz 

oo 



= E l^t^^"""^'^ + r^^'S'{r))S\r){r + p)-('+^)] . 

Formula fllO.lip shows that in order to compute the contribution to the density of 
states from H(Q3)p, we need to integrate the RHS of (110. ISp against dX (or rather du) 
and We are going to integrate against d^ first. We will prove that this integral is a 
convergent series of products of powers of p and In p. The coefficients in front of all terms 
will be bounded functions of X, so afterwards we will just integrate these coefficients to 
obtain the desired asymptotic expansion. 

Let us discuss how the RHS of (110. ISp depends on the coordinates X and $ (or rather 
$). Equations flOTTOD - flOTD . Lemma [73] and Remark [TO show that the RHS of f[To:T8|l 
is a sum of terms of the following form: 

(10.19) C///i(X)/2($)/3(X;p;$). 

Here, f\ is a uniformly bounded function of X coordinates only. It consists of contri- 
butions from the cut-off functions (fe with G 2J and from the terms in (I9.20p . (I9.2ip 
corresponding to OpO'j e 03. The function /2($) is a product of powers of {sin$q}. 
This function comes from differentiating (17.151) with respect to r. Finally, /a is of the 
following form: 

T 

(10.20) MX; p; $) = l[{k + pJ^K M%)r''- 

t=l g 

This function corresponds to the negative powers of inner products dt) given by 
Lemma [731 part (ii). Here, {6*} are coefficients in the decomposition {6t)^± = J^q^qP'q'^ 
recall that these numbers are all of the same sign and satisfy (I7.16p . Without loss of 
generality we will assume that all 6* are non-negative. The number It = l{b\, . . . , '■— 
{X,{Ot)w) + Lm+i^qh\ satisfies Pn""^^pV < < Pn'^^Vn'^; siucc our assumptions 
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imply \{X,0(^)\ ^ p"™. This number depends on X, but not on $ or p. The number 
kt = k{b\, . . . , b^K+i) is positive, integer, and independent of ^. Our next objective is to 
compute the integrals of (110. 19p over the domain {$ G Mp} and prove that these integrals 
enjoy asymptotic behaviour (13. 9p with uniformly bounded coefficients (as functions of X). 
The calculations will be rather messy technically, although the main ideas of computing 
them are not too difficult. 



10.2. Computing the model integral. Before computing the integral of (110. 19p . we 
will deal with a simpler integral 

(10.21) Jk - ' ' ' ^1 •••^i^ a^?,...a^?K 



^0 



and then we will discuss how to reduce our initial integral to (110. 2ip . Here, nj, kt, k[ G 
(N U {0}), 7 < 1, < P"^+°~ < < P"''+°+ < fhl\ < Q < f^r!\ p~^" < 6* < pfl\ 
and p^^" ^ 6* ^ p^, where 5o > is sufficiently small (for the sake of definiteness, 
we put 5q := grfgQQQ ; obviously, we assume that these inequalities hold only for non-zero 
values of 6* and 6*). We introduce the following notation: 

(10.22) i^:=E^^' Q--=J2^^^ Q'--=Y.^'t 

j t t 

and we sometimes will denote the integral (110. 2ip as Jk{P, Q, Q')- We will also need the 
auxiliary positive numbers pj, g^, j = 0, . . . , (i, defined by 

^^=3^00' ^^ = ^^ + 340' ^-^=^^+^^ + 34oO='^- 
Obviously, po < 1/100. 

Lemma 10.2. Assume that Ct ^ Pn'^^- Then we have: 

K oo 

(10.23) Jk{p, q, q') = Y.^Hp)r J2 Q')p~'^ 

q=0 p=0 

where 

T 

(10.24) \e{p,q-P,Q,Q')\ « p(f/3-^-)^p;«^2«' 

t=i 

These estimates are uniform in the following regions of variables: 



< 7 < 1, p(; < < p„^ , p;'"-' < Q < p 



n 5 



'10 25) 

^ ■ p;^" « 6* « p^ p;^" « 6* « p^ p^f < p. 
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Remark 10.3. The estimates (110.241) are more natural than they may look. Indeed, 
each time we run the inductive argument in the proof (i.e. each time we increase K), we 
apply the geometric series expansion, which results in a slight worsening of the estimates. 
This accounts for the need to have px in the exponent of p„. 

Proof. The proof will go by induction in K. The base of induction [K = 0) is trivial 
(and the case K = 1 has been discussed in pi]). Suppose, we have proved this statement 
for i^T = 5* — 1, and let us prove it for K = S. 

Step I. First, we consider the area where $5 > p^^"^. We do not change terms in 
the denominator of Jk where 6^ = 0. If 65 7^ then we proceed with the following 
transformations : 

(10.26) 

where constants Cm{kt) satisfy the estimate 

(10.27) \CM\ < pTI"^^^' " ^1 < P^'2'"+'='-i. 

\ m J 

Now, we can move powers of p out of the integral and denote q := h^g^s- Obviously, q 
satisfies (110. 25p with K = S — \. For terms which do not contain p we just denote q := 
Ct + h^s^s- Then, we can apply the induction assumption for K = S — 1. Corresponding 
coefficients will depend on $5 uniformly. As a result, we obtain the following expression 
as the contribution to Jk from the region {$5 > p^^^} (we denote m := mi + ■ ■ ■ + uit 
and assume for simplicity that 7^ for all t): 
00 00 

(10.28) E E C^^{ki)...C^,{kT)p-'^-'^Js-i{P.^.m + Q + Q'). 

mi=0 mT=0 

Of course, each time we write Js-i{P,0,m + Q + Q'), it denotes a different integral of 
the form (110.211) . but by the assumption of induction all of them satisfy 

5-1 00 

(10.29) Js-i{P, 0,Q + m + Q') = Y,iHp))' E ^^P' 0, Q + m + Q')p-P 
with 

T 

(10.30) \e{p, g; P, 0, Q + m + g')l < p(^/^-Ps-i)PpiQ+^)'io2Q+^+Q' J] c^^* 

t=i 
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(when we use f ll0.24p as the induction hypothesis, we replace Ct and q by the corre- 
sponding lower bounds p~'^° and q). 

Contribution from this region of integration into coefficients e{p, q) of the integral Js 
can therefore be estimated from above by the following expression: 

(10.31) 

P-Q ( T \ 

^(2/3-ps-i)(p-m-Q)^(Q+m)<?o^m/22m+Q-T2Q+»n+Q' | "["[" c^^* ) ^ 1< 

m=0 \t=l / ra\-\ h»nT=m, mj>0 



hi 



X 



Q/3+Q?0~(2/3-PS-l)Q22Q-^ ^ ^mgo+m/2-(2/3-ps-i)m22™2™+'^"-^ < 



m=0 



^(2/3-..).^-c5/.2Q'JJe7^^ 

Notice that at this step we have S" — 1 as the largest power of ln(p). 

Step II. From now we are in the area $5 < p^^^ ■ Then we can transform all terms 
not containing p in the denominator of Jk'- 

(10.32) (^ 65$, + cY^'^ = ic,)-''^ (l + ' = CiiK) 



Ct 



m=0 \ j 



where 

(10.33) \C'„XK)\ < c7^'""'2'"+^'-^ < pf'^c7^*2'"+'='-\ 

Then, only terms with p are left in the denominator. We change variables xj := ^jp and 

obtain the integral 

(10.34) 



First, we consider the integral along the region where < xs < fh!^ '^^ ^ Obviously, 
the corresponding contribution to Js can be computed as 



00 

p-P-sj-p-C. 

m=0 

where 



|L^m| ^ Pn Pn ^ Pn ^ Pn X 1 ' 



T 

it' 
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Thus, if we put m = p — P — Swe obtain the following estimate for the coefficient in 
front of (notice that ps < qs-i — Qs ^ 25o for 5* > 1): 

T 

(2/3-ps)p -Q/3 -(P+5)gsoQ' TT ^-^^ 
Pn Pn Pn ^ J_ J_ • 

t=l 

Step III. The case when xs € (pn^^ ''^"^pp"^^). Once again if 6^ = then we leave 
such terms unchanged. If 65. 7^ and thus 65 > p^''°, we perform the following transform 

to ffTo^D (cf. (daMD): 



S 00 / ^ 

(10.35) {k + J2b]^jr'' = Y.Cm{k,) — - 



m+kt 



j=i m=o v^i^i-^J, 

and introduce new variables Zj := Xj/xs, j = I, . . . , S — 1. Thus, we reduce the problem 
to the integrals of the following form: 



Now we can remove xs from the second bracket in the denominator and then apply 
the induction assumption for the internal 5 — 1 integrals (i.e. the integrals against 
dzi . . . dzs-i with q := 65 and p := xs- This induction assumption guarantees that 
these internal integrals can be expressed as a series in powers of a;^ and Inxs, with the 
biggest power of Inxs' being 5 — 1. Then, we multiply this expansion by a (possibly 
negative) power of xs and integrate the product against dxs- As a result, we obtain 
a decomposition (110.231) of Js, with the biggest power of Inp being equal to S. The 
estimate of the contribution of Step III to the coefficients e{p, q) is similar (but rather 
more tedious) to the estimates in the first two steps, and we will skip it. □ 

10.3. Reduction to the model integral. Now we will discuss how to deal with our 
initial integral 



(10.36) Jk:= J 



(sin$i)"i . . . (sin$/^)"'ff(sin <I>k+i)"^+^ rf$ 

ul,ih+pj:f=ibUM,r^ • 



The main problem with reducing the integral along Mp (or even along Mp fl {$1 < 
• • • < $^ < ^K+i}) to the model integral (110. 2ip is the limits of integration: the upper 
limit of integration against d^x is not a constant (since the collection of points where 
$A' = ^K+i has variable coordinate $/<). In order to rectify this, we define 

(10.37) M := {sin $1 < ■ ■ • < sin <l>,_i < p-P'*, p-f'* < sin q = s,...,K+ 1}. 

It is clear that Mp can be represented as a union of several domains of this type. Lemma 
17.41 shows that we always have at least one 'large' variable in M, i.e. s < K +1. We also 
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introduce the 'spherical' coordinates in the . . . , $A'+i)-subspace: we put 

sin% = <lj-, j = l,...,s-l, 

sin $s — p^^'' = f cos <l>s, 

sin $5+1 — = f sin <l>s cos 

(10.38) ^ ^ ' 

sin — Pn^"^ = r sin sin $^+1 ... sin ^k-i cos 
sin ^K+i — Pn^'' = r sin sin <E>s+i ... sin ^k-i sin <l>x, 

so that (f, <l>i, . . . , are the new coordinates. Since only K of the coordinates $j were 
independent, we can consider the new variables ($i, . . . , as independent (and f as 
a function of the independent variables). We remind that (see Lemma [7.41) we always 
have maxj sin $j ^ and thus 

(10.39) < f < 1. 

The point in introducing these variables is that the limits of integration over M become 
simple: 

(10.40) / d^K d^K-i d^J d^s-1 d^s-2.-- d^,. 
Jo Jo Jo Jo Jo Jo 

When we insert the values of sin $g, q = 1, . . . , + 1 given by (110. 38p into (I7.12p . we 
obtain the quadratic equation for finding f: 

(10.41) i,f2 + 2Bsf + {Cs - 1) = 0, 
where 

As = '^{ajs cos + ajs+i sin cos "l^+i + . . . 



10.42) 



(10.43) 



and 



3 

+ aj K+i sin <l>s sin ^s+i ... sin ^k-i sin ^x)^ > 0, 

s-l K+1 

Bs = (^jq^q + Pn'''' Cijq){<^js COS + ttj s+1 siu $5 COS ^ s+l + • • • 

j q=l q=s 

+ ttjK+i sin sin <l>s+i ... sin <lx-i sin <lx), 



s-l K+l 

:io.44) = Y^iY. ''n^'^ + pn"' E «^'?)' > 0- 

3 <?=! 9=* 
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Therefore, we have 



-B,±xlB^,-AsCs + A, 
(10.45) r = 



As 

Note that Cauchy-Schwarz inequahty imphes AsCs > B^. Since equation (I10.4ip obvi- 
ously has at least one real solution, we have < As + B^ — AsC s < As. Next, due to 
Lemma [7.31 we have 

(10.46) \Bs\ « \Cs\ « p-'"''+°+. 
Then using (IT(09l) . ( TTOiTll . and (ITIM we get 

(10.47) pr » ^ > > ^ » 1- 

Since f is positive we obviously have 



Bs + - AsCs + A, -S, + ^Asdl - (C, - 



and thus f is analytic with respect to -B^, (7^, i.e. with respect to all $j, j = 1, . . . , s — 1, 
uniformly in / = s, . . . , i^', inside M. It is easy to see from fl7.13p . f l7.14p and f llO.SSp 
that the same is true for the Jacobian 



We also notice that, if we denote 
(10.48) fo := r - ij^/^ 

then fo satisfies the same analyticity properties as f and fo = 0(p„^''+°^). 
Thus, we arrive at the integrals of the following form: 

7r/2 p-K/2 pTT/2 



'10 49) 



d^s-i / d^s-2--- / d^i- rj. ... 

Here the function -F($s, . . . , ^k) is uniformly bounded with respect to <l>s, . . . , in M 
and 

s-l iS'+l 

(10.50) S = ^($1, ...,$;,):= J] 65$,. + ^ fej.p-^'* + (i;V2 + fo)F($„ . . . , 1-^), 

i=i j=s 

where 

F := 6* cos + sin cos l>^+i H h sin $5 sin $3+1 ... sin l>x-i sin l>x. 
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1 



foF 



Now we apply the construction from the STEP I of the proof of Lemma 110.21 We do 
not change terms {k + pS) with 6* = for all j = s, . . . , i^' + 1 (such terms are equal to 

{k + P^jllb^j^j)). Otherwise, we write (cf. f lTCT]) ) 

m+kt 

ik + pS)'"^ = J2CUkt 

It remains to notice that 

s-l K+l / 

s~' = iyb'A+yb]p-'''+A:'/^FrUi + —^ . . 

h k \ ErJi^'^^+Ei^'b^^Pn'^+Aj'/'F^ 

and decompose the last expression using geometric progression. We remind (see fll0.47p 
and (110. 48p ) that foAl^"^ <C p~^''~^^~^. Since fo is analytic in $i, . . . , ^^-i, we end up with 
the model integrals Jg-i (with q := Yl!j=s b^jPn^"^ + Ag^^'^F) uniformly depending on 
the parameters <l>s, . . . ^k+i- Summing this and using Lemma [10. 2[ we have proved the 
following result. 

Lemma 10.4. We have: 

K oo 

(10.51) Jk = J](ln(p))^J]e(p,g)p-^ 

q=0 p=0 

where 

(10.52) |e(p,g)|«pi'/=^-^-)V-«'^. 

These estimates are uniform in the following regions of variables: 

(10.53) p^«/i«py2^ p;^°«6;.«p^°, p^/'-^-<p. 

Now Lemma 110. 4^ Remark 110. H and equation f llO.lSp show that the integral (110. lip 
admits decomposition of the form (110. 5ip for 0.7p„ < p < 5pn- This, together with 
equations (110. lip . (110. 2p . (I6.16p . Corollary 14.31 and the observation that the number of 
different quasi-lattice subspaces QJ is < p^"*", completes the proof of Lemma I3.3l and. thus, 
of our main theorem in the case of all domains Mp being simplexes. It remains to discuss 
how to reduce the case of general region S(QJ)p to the case of a simplex. 

11. Integration in non-simplex domains 

Now let us consider the case when the number Jp of defining hyperplanes is bigger 
than K + l. Recall that we have 

(11.1) 3(QJ)p = G G fi(QJ) & {$^^,^ij{p)) > W, J = 1, • • • , Jp}- 

In this section, we give only a brief description of how to prove the main statements, 
since the complete proof would be too long and tedious. The complexity of the proof is 
mostly caused by the fact that we need to describe a procedure working in all dimensions. 
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If one is interested only in the cases d = 2 or d = 3, the proofs become substantially 
simpler (indeed, the case d = 2 has already been proved, since then we have i^' = 1 or 
K = 0, and any polyhedron in dimensions or 1 is obviously a simplex). 

Consider first the case when S(5J)p is a 'cone', i.e. there is a point a = a{p) G ^J-*- 
such that (a, fijip)) = Lm+i, j = I, . . . , Jp (such point, if it exists, is always unique, and 
the existence of it is automatic in the simplex case, i.e. when Jp = K + 1). Then we can 
introduce the coordinates (r, $) by formulas (17.61) and (17.71) with $ G Mp, where Mp is 
still given by (17. 8p . 

Lemma 11.1. Suppose, G 0^,. Then we can write 6<xii^ = Y2q^qi^qip) > where either 

all hq are non-positive, or all of them are non-negative (hut such a decomposition is not 
necessarily unique). 

Proof. Our assumptions imply that for each r] G S(5J)p (where S(QJ)p is defined by 
(17.41) ) we have {r],6y±) ^ 0. Assume for definiteness that {r},Oy±) > 0. This property 
can be reformulated like this: whenever z G 2?"*" is a vector with {z,6^±) < 0, there is 
at least one vector /ij(p) such that {z,fij{p)) < 0. This, in turn, is equivalent to saying 
that whenever z G QJ"*" is a vector with (z,0(jj±) > 0, there is at least one vector fj,j{p) 
such that {z,fj,j{p)) > 0. Consider the set S := {z = J2q^qP'qiP)'^q — 0}- need 
to prove that 0<jj^ ^ S. Suppose not. Let zq be a nearest to 6<^± point from S. Then 
{0^± — zo,0Qjx) > 0, because otherwise the point \0^±\'^ {zq, 6<^±)~^zo G S" is closer to 
d<^± than Zq. Thus, there is a value of j, say j = 1, so that (^sgx — zq, f^iip)) > 0. But 
then for sufficiently small e > the vector zq + efj,i{p) G S" is closer to 0qj± than zq. This 
contradiction proves our lemma. □ 

This result shows that there is a big similarity between cases of the cone and the 
simplex. The only difference from the simplex case is that the number of sides of Mp is 
now greater than K-\-l. This however means that if we introduce the angular coordinates 
$ in the same way as in Section 7, we will have difficulties trying to get rid of several of 
them. Instead, we will follow a different strategy: we will cut the spherical polyhedron 
Mp into several simplexes: 

(11.2) Mp = Uq^'Mp 

and then perform integration over each simplex '^Mp in the same way as we did in sections 
7-10. The only thing we need to make sure is that the lengths of all sides (edges) of '^Mp, 
as well as all non-zero angles between two sides of any dimensions of "^Mp, is ^ (let 
us call this angles and sides property). We, however, know what that angles and sides 
property holds for the original polyhedron Mp because of Lemma 17.11 (strictly speaking. 
Lemma [7. II was proved for simplexes, rather than for cones, but the proof is the same). 
Thus, the only problem we face is how to cut a polyhedron Mp into simplexes '^Mp 
without drastically decreasing sides or angles. We do it by induction in K. For K = 1 
the statement is obvious (each 1-dimensional connected polyhedron is a simplex, i.e. an 
interval). Assume that K is arbitrary. Also assume for simplicity that Mp is not a 
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spherical, but a Euclidean polyhedron (we can achieve this by projecting Mp onto any 
hyperplane tangent to it; obviously, this projection keeps the angles and sides property 
invariant). 

Step I. We find a simplex Mp C Mp satisfying the angles and sides property. To do 
this, we consider a ball centered at any vertex v of Mp of radius ^ but sufficiently 
small so that the intersection of this ball with Mp is a cone. The intersection of the 
boundary of this ball with Mp is a polyhedron Np of dimension K — 1. Running the 
induction argument, we can find a {K — l)-dimensional simplex Np C Np satisfying the 
angles and sides property. Now we define Mp as a convex hull of v and vertexes of Np. 
A straightforward geometrical argument implies that Mp satisfies the angles and sides 
property. In particular, the volume of Mp is ^ 

Step II. We find a point 77* G Mp such that the distance from r)* to each of the 
{K — l)-dimensional sides of Mp is ^ (the distance to the side is the length of the 
perpendicular dropped to the hyperplane containing this side). This point can be chosen 
to be the center of gravity of Mp. Indeed, the distance from 77* to a {K — l)-dimensional 
side of Mp is the average of the distances from the vertexes of Mp to this side. Thus, 
we need to show that the distance from at least one of the vertexes of Mp to this side is 
^ But if this were not the case, then the breadth of Mp in the direction orthogonal 
to that side were ^ p°~, and so the volume of Mp were <^ p°~, which would contradict 
estimates from Step I. 

Step III. Now we use the inductive assumption and cut each {K — l)-dimensional side 
of Mp into simplexes. Taking convex hulls of 77* with these simplexes, we obtain the 
required decomposition of Mp into simplexes '^Mp. It is a geometric exercise to check 
that such constructed simplexes ''Mp satisfy the angles and sides property. 

Let us denote by "^3(51)^ the infinite cone with the vertex a and a cross-section ''Mp, 
i.e. 

(11.3) ''H(QJ)p := a + {I G G fi(5J) & Yv{i^^) G "Mp}. 
Then we obviously have 

(11.4) S(QJ)p = U,''H(2J)p. 

Now let us discuss how to perform the integration over ''S(QJ)p. Let us fix q and p and 
denote by i>i, . . . ,i>k+i the interior unit normal vectors to the faces of ''3(2J)p. We 
denote, as before, $q := | — 0(^fr;± — a, Vq{p)), q = 1, . . . , K + 1. In order to perform 
the integration, we need to check that Lemma 17.51 is still valid in the cone case. So, let 
6 G 0^. Applying Lemma flLlj we deduce that 

(11.5) d^. = Y,KKip)^ 

q 

where either all bg are non-positive, or all of them are non-negative; assume for defi- 
niteness that all of them are non-negative. We also have (applying, for example, the 
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same lemma) that each vector fiq{p) admits a decomposition fiq{p) = J2h=i^ Ki^i with 
all coefficients bqi being non- negative. Now we have (denoting, as usual, rj := $,<^± — a 
and putting k := ^^L, bqbqi > 0, / = 1, . . . , + 1): 

0) = (X, 0^) + (a, 02jx) + (ry, 02jx) 

Jp K+l 

Q^ = far) + Yl ^"(a, A,) + 5^ kirj, i^i) 

^ ' ^ q=l 1=1 

= {X, 6>5j) + IqLm+i + r ^ 6; sin 

The estimates < \bq\ < p^^ can be proved in the same way as Lemma [2.81 Finally, 
fill. 51) implies that Ylq K ^ Pn^ 2J. Multiplying flll.Sp by a, we deduce that 

This finishes the proof for the cone case. Now let us discuss the general case. 
Let a be any point inside 3(2J)p such that for all j we have L^+i < {sl, flj{p)) ^ 
Lm+iP*^- For each / = 0, . . . , Jp, we define 

(11.7) (£2Ti,Aj-(p)) > (a,/i^-(p)), J = 1,...,/, & 

(^2Ji, Aj(p)) > ^m+l, j = / + 1, . . . , Jp} 

and 

(11.8) s{^)i:=Emi\B{^y;\ 

Then we obviously have 

(11.9) H(^)p = Ei^Y/ U (uJo ^H(QJ);,) 

(as usual, modulo boundary points). The domain H(QJ)p'' is a cone, so we already know 
how to deal with it. Now let us consider 3(QJ)p for some /. Let us introduce a coordinate 
t = t{$,) := (^, /i;_|_i(p)). Then for ^ G S(2J)p we have t G [Lm+i, (a, /i;_|_^(p))]. For each 
t G [L^+i, (a,/ii+i(p))] we denote 0\t) := {$, G H(2J)|,, (^,/i;+i(p)) = t}. It is easy to 
see that the domain 0'(t) is of the same type as the domain S(5J)p, i.e. it can be written 
as 

(11.10) {^gM^ ^^Gl](5J),(|,Am(p))=^ & + 1 {^^^,i>j)>s,}, 

where Uj are, of course, normalized projections of some /ij(p) onto the plane orthogonal 
to iiij^^{p). Our aim is to compute 

(11.11) voli+nH(Q3)'-voli"nH(QJ)' 
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(or at least prove that this expression admits an asymptotic expansion in powers of p 
and In p) . But flll.lip is obviously equal to 

(11.12) / {yo\{a+ r\0\t))-Yo\{A- no\t)))dt. 

J Lm+\ 

If dimO'(t) = 1, or, more generally, if 0\t) is a simplex, then we can perform integration 
over 0\t) as described above, since the formula (I11.6P would still be valid, in the sense 
that 

(11.13) (1,0) =C(X,t)+r^6zsin$z, 

I 

where Lm+ipt' < C{X,t) < Lm^ip^ and the coordinates (r, {$;}) are the shifted polar 
coordinates in 0\t). Results of Section 10 show that for each fixed t the expression 
(vol(A+ n 0\t)) — vol(A~ n 0'(t))) admits the asymptotic expansion in powers of p 
and Inp, with coefficients being uniformly bounded in t (and X). Now it remains to 
integrate this expansion against dt (and dX). If 0\t) is a cone, we cut it onto simplexes 
as described earlier in this section, and then integrate over each simplex separately. 
Finally, if 0\t) is not a cone, we continue the process of reducing dimension until the 
dimension of 0\t) becomes equal one. 



References 

G. Barbatis and L. Parnovski Bethe-Sommerfeld conjecture for pseudo-differential perturbation, 
Comm.P.D.E., 34(4) (2009), 383-418. 

B. Helffer, A. Mohamed, Asymptotics of the density of states for the Schrodinger operator with 
periodic electric potential, Duke Math. J. 92 (1998), 1-60. 

M. Hitrik, I. Polterovich, Regularized traces and Taylor expansions for the heat semigroup, J. London 
Math. See. 68(2) (2003), 402-418. 

M. Hitrik, I. Polterovich, Resolvent expansions and trace regularizations for Schrodinger opera- 
tors. Advances in Differential Equations and Mathematical Physics, Contemporary Mathematics, 
American Mathematical Society, 2003. 

Yu. Karpeshina, On the density of states for the periodic Schrodinger operator, Ark. Mat. 38 (2000), 
111-137. 

Yu. Karpeshina, Perturbation theory for the Schrodinger operator with a periodic potential, Lecture 

Notes in Math., Vol. 1663, Springer Berlin 1997. 

T. Kato, Perturbation Theory for Linear Operators, Springer 1980. 

E. Korotyaev and A. Pushnitski, On the High-Energy Asymptotics of the Integrated Density of 
States, Buh. LMS 35 (2003), No. 6, 770-776. 

M. A. Naimark Normed rings, P.Noordhoff N.V., Netherlands, 1959. 
L. Parnovski Bethe-Sommerfeld conjecture Annales H. Poincare, 9(3) (2008), 457-508. 
L. Parnovski, R. Shterenberg Asymptotic expansion of the integrated density of states of a two- 
dimensional periodic Schroedinger operator, Inv.Math., 176(2) (2009), 275-323. 
L. Parnovski, A. Sobolev Bethe-Sommerfeld conjecture for periodic operators with strong perturba- 
tions, arXiv:0907.0887vl, submitted for publication. 

M. Reed, B. Simon, Methods of Modern Mathematical Physics, Academic Press, 1978. 



54 



L. PARNOVSKI & R. SHTERENBERG 



[14] A. V. Savin, Asymptotic expansion of the density of states for one- dimensional Schrddinger and 
Dirac operators with almost periodic and random potentials, Sb. Nauchn. Tr. , I.F.T.P., Moscow 
(Russian), 1988. 

[15] M. Shubin Almost periodic functions and partial differential operators, Russian Math. Surveys 33(2) 
(1978), 1-52. 

[16] M. Shubin The spectral theory and the index of elliptic operators with almost periodic coefficients, 
Russian Math. Surveys 34(2) (1979), 109-157. 

[17] D. Shcnk and M. Shubin, Asymptotic expansion of the state density and the spectral function of a 
Hill operator, Math. USSR Sbornik 56 (1987), No. 2, 473-490. 

[18] M. Skriganov, Geometrical and arithmetical methods in the spectral theory of the multi- dimensional 
periodic operators, Proc. Steklov Math. Inst., Vol. 171, 1984. 

[19] A. V. Sobolev, Asymptotics of the integrated density of states for periodic elliptic pseudo- differential 
operators in dimension one. Rev. Mat. Iberoam. 22 (2006), no. 1, 55-92. 

[20] A. V. Sobolev, Integrated Density of States for the Periodic Schrddinger Operator in Dimension 
Two, Ann. Henri Poincare 6 (2005), 31-84. 

[21] O. A. Veliev, On the spectrum of multidimensional periodic operators. Theory of Functions, func- 
tional analysis and their applications, Kharkov University, 49 (1988), 17-34 (in Russian). 

Department of Mathematics, University College London, Gower Street, London, 
WCIE 6BT, UK 

E-mail address: LeonidQmath.ucl.ac.uk 

Department of Mathematics, University of Alabama at Birmingham, 1300 University 
Blvd., Birmingham AL 35294, USA 
E-mail address: shterenb@iiiath.uab.edu 



